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Abstract. In this review I discuss the application of scanning magnetic imaging to 
fundamental studies of superconductors, concentrating on three scanning magnetic 
microscopies - scanning SQUID microscopy (SSM), scanning Hall bar microscopy 
(SHM), and magnetic force microscopy (MFM). I briefly discuss the history, sensitivity, 
spatial resolution, invasiveness, and potential future developments of each technique. I 
then discuss a selection of applications of these microscopies. I start with static imaging 
of magnetic flux: An SSM study provides deeper understanding of vortex trapping in 
narrow strips, which are used to reduce noise in superconducting circuitry. Studies 
of vortex trapping in wire lattices, clusters, and arrays of rings and nanoholes show 
fascinating ordering effects. The cuprate high-T c superconductors are shown to have 
predominantly d-wave pairing symmetry by magnetic imaging of the half-integer flux 
quantum effect. Arrays of superconducting rings act as a physical analog for the Ising 
spin model, with the half-integer flux quantum effect helping to eliminate one source of 
disorder in antiferromagnetic arrangements of the ring moments. Tests of the interlayer 
tunneling model show that the condensation energy available from this mechanism can 
not account for the high critical temperatures observed in the cuprates. The strong 
divergence in the magnetic fields of Pearl vortices allows them to be imaged using 
SSM, even for penetration depths of a millimeter. Unusual vortex arrangements occur 
in samples comparable in size to the coherence length. Spontaneous magnetization 
is not observed in Sr2Ru04, which is believed to have p x ± ip y pairing symmetry, 
although effects hundreds of times bigger than the sensitivity limits had been predicted. 
However, unusual flux trapping is observed in this superconductor. Finally, unusual 
flux arrangements are also observed in magnetic superconductors. I then turn to 
vortex dynamics: Imaging of vortices in rings of highly underdoped cuprates places 
limits on spin-charge separation in these materials. Studies of spontaneous generation 
of fluxoids upon cooling rings through the superconducting transition provide clues to 
dynamical processes relevant to the early development of the universe, while studies 
of vortex motion in cuprate grain boundaries allow the measurement of current- 
voltage characteristics at the femtovolt scale for these technologically important 
defects. Scanning SQUID susceptometry allows the measurement of superconducting 
fluctuations on samples comparable in size to the coherence length, reveal stripes in 
susceptibility believed to be associated with enhanced superfluid density on the twin 
boundaries in the pnictide superconductor Co-doped Ba-122, and indicate the presence 
of spin-like excitations, which may be a source of noise in superconducting devices, in 
a wide variety of materials. Scanning magnetic microscopies allow the absolute value 
of penetration depths to be measured locally over a wide temperature range, providing 
clues to the symmetry of the order parameter in unconventional superconductors. 
Finally, MFM tips can be used to manipulate vortices, providing information on flux 
trapping in superconductors. 
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1. Introduction 

Superconductivity, the complete absence of electrical resistance in some metals below a 
critical temperature T c , is one of the best understood phenomena in solid-state physics. 
As described in the Bardeen-Cooper-Schrieffer (BCS) theory [T], superconductivity 
occurs in conventional superconductors through phonon-mediated pair-wise interaction 
of charge carriers to form Cooper pairs. The characteristic energy needed to form a 
single particle excitation from the superconducting condensate is the energy gap A. 
The Cooper pairs form a macroscopic quantum state characterized by a wave function 
ip{r) = |"0(r)|e*^^, where <fi is the Cooper pair phase, with the density of Cooper pairs 
n s (r) = ^(r)] 2 . The macrosopic quantum nature of the pairing state leads to a wide 
variety of behaviors: persistent currents, Meissner screening, flux quantization, the 
Josephson effects, and superconducting quantum interference. Superconductivity has 
become a multi-billion dollar industry, with applications in medical imaging, cell phone 
filters, electrical motors, power transmission, transportation, magnetic manipulation, 
magnetic sensing, etc. 

However, there is still much to learn about superconductivity: Although 
superconductivity in the copper-oxide based perovskites [2j [3] was discovered in 1986, 
despite almost 25 years of intense effort a consensus on the pairing mechanism in the 
cuprates has yet to emerge. The heavy Fermion superconductors [UE] have large carrier 
masses, strong interaction between spin and charge degrees of freedom, and potentially a 
wide variety of Cooper pairing symmetries. The non-cuprate perovskite superconductor 
Sr2Ru04 [6] is believed to have a p-wave pairing state that breaks time reversal symmetry 
[7] . Interest in unconventional superconductors has been reawakened with the discovery 
of iron-based compounds with high critical temperatures [8]. 

Scanning magnetic microscopies have played, and will continue to play, an 
important part in our efforts to understand superconductivity. There are now a large 
number of scanning magnetic microscopies. Although it is beyond the scope of this 
review to cover all of these techniques in detail, it is appropriate to briefly describe 
some, and provide a few citations as an introduction to the literature: 

In Lorentz microscopy |H] QUI EEH H2], a coherent beam of high energy (300keV) 
electrons are passed through a thin sample and then imaged slightly off the focal plane. 
Magnetic fields introduce small phase shifts which cause interference in the defocussed 
image. Lorentz microscopy can image vortices with deep sub-micron spatial resolution 
at video rates and is sensitive not only to surface fields but also to fields within the 
sample. However, the sample must be thinned. 

In magneto-optic microscopy [I3J [HJ [T5], the sample is illuminated with linearly 
polarized light. It passes through a magneto-optically active layer, reflects off the sample 
surface, passes again through the magneto-optically active layer, and then through a 
crossed polarizer. Regions of the sample with magnetic fields produce rotations of 
the polarization axis of the light, which is detected as a bright spot. Magneto-optic 
microscopy has the advantage of relative simplicity, but has limited spatial resolution 
and sensitivity. Nevertheless, observation of individual superconducting vortices with 
0.8/xm spatial resolution has recently been reported [T5] . 

In Sagnac interferometry [161 EL71 EEE], a single beam of light is split into two 
components which travel along identical paths in opposite directions around a loop. 
Any effects which break the time reversal symmetry, including the magneto-optical 
effect, will cause interference, which can be sensitively detected. High spatial resolution 
can be achieved using scanning near-field optical microscopy [T7] . 

In spin polarized scanning tunneling microscopy [TjJl EH EI] a scanning tunneling 
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microscope has a tunneling tip with a ferromagnetic or antiferromagnetic atom from 
which the tunneling electrons emerge. The tunneling probability depends on the 
spin polarization of the sample, and therefore contrast is achieved if there are 
spatial variations in the atomic spin polarizations. Spin polarized tunneling is 
useful for studying superconductivity in, for example, the high-temperature cuprate 
superconductors, since it is believed that magnetism and superconductivity are 
intimately related in these materials. However, the spin polarized tunneling current 
is not directly related to the local magnetic fields. 

In diamond nitrogen- vacancy microscopy [22j [23j EH EHJ [26], a magnetic field 
dependent shift in the energy levels of a single nitrogen-vacancy center in diamond 
is detected optically. There are two basic schemes for detecting this energy shift. In 
the "dc" technique, microwaves at a fixed frequency as well as light are incident on 
the center. The intensity of the resultant fluorescence depends on the local magnetic 
field. This method can detect milliTesla fields with, in principle, nanometer scale spatial 
resolution. In the "ac" technique, spin-echo microwave techniques are used to detect 
the modulation of the center's energy levels by the local magnetic field. The "ac" 
technique has a sensitivity of 60 nT, with in principle a spatial resolution of a few nm, 
but has reduced sensitivity at zero frequency. These techniques have the advantages of 
good sensitivity, small invasiveness, room temperature operation, inherently high speed, 
and potentially high spatial resolution. It remains to be seen whether single nitrogen- 
vacancies in nm sized diamond crystallites can be placed at the end of scanning tips 
while still retaining long lifetimes. This is necessary to achieve high spatial resolution 
and sensitivity simultaneously. 

A magnetic tunnel junction [271 123 [29] is a planar structure with two ferromagnetic 
electrodes. The tunneling current through the junction depends on the relative 
alignment of the magnetic moments of the electrodes, and therefore depends on 
the magnetic field environment. Magnetic tunnel junctions suitable for scanning 
microscopy are commercially available. They currently have lateral dimensions (and 
therefore spatial resolutions) of about 4 microns, and field noises at 100 Hz of about 
5x 10~ 8 T/Hz 1 / 2 [2H]- They have the advantage of room temperature operation, but are 
not currently as sensitive as Hall bars or SQUIDs. 

There are already excellent reviews covering the topic of magnetic imaging of 
superconductors [3U1 EE E2], and there has been an enormous amount of work in 
this area. I will therefore not attempt to survey the field exhaustively, but will only 
cover recent work in three scanning magnetic microscopy techniques: scanning SQUID, 
scanning Hall bar, and magnetic force microscopies. I will present a short introduction 
to the basics, advantages and shortcomings of each technique, and then present a 
selection of fundamental applications of these magnetic microscopies in the area of 
superconductivity. 

2. Fundamentals 

2.1. Scanning SQUID microscopy 

Josephson [33J predicted in 1962 the possibility of coherent tunneling of Cooper 
pairs between two weakly coupled superconductors, and wrote down the fundamental 
relations: 



I s = I sin ip 

h dtp 



where I s is the supercurrent through the weak link, ip is the difference between the 
Cooper pair phases on the two sides of the weak link, V is the voltage across the weak 
link, and Jo is the critical current, the maximum current that can pass through the 
weak link before a voltage develops. A physical Josephson junction can be modeled as 
an ideal Josephson element in parallel with a capacitor and a resistor in the resistively 
shunted junction (RSJ) model [MIES]- This model obeys the equation: 

h dip hC d 2 p> 

lB = ^RlI + IoSm ^ + ^^ + In > (2) 
where Ib is the bias current through the weak link, R is the resistance, C is the 
capacitance, and I n is a noise current, given for the case of an ideal Johnson noise 
source resistor by 

<I n > 2 =4k B TAf/R, (3) 

where < I n > 2 is the time-averaged current noise squared, and Af is the frequency 
band width over which the current noise is measured. This equation has the form 
of the equation of motion of a driven, damped harmonic oscillator. The dynamics of 
a Josephson weak-link can be visualized as that of a particle moving in a sinusoidal 
"washboard" potential with oscillation amplitude l <& /27r, where $o — h/2e is the 
superconducting flux quantum, and the average slope of the potential is given by 
< dE/dip >= $o(Js — Jn)/27r. In the absence of a noise current, the particle will 
escape from a local potential minimum and run down the washboard potential if the 
bias current Ib is greater than I , in the process developing a voltage given by the second 
of Eq.s[TJ The solutions to Eq. [2] can be divided into two classes, according to whether 
the Stewart-McCumber parameter /3 C = 2ttIqR 2 C/^o is greater or less than 1. If /3 C < 1, 
a running particle retraps into a local potential minimum (V = 0), as soon as Ib = Iq, 
and the current- voltage characteristic is non-hysteretic. On the other hand, if f3 c > 1, a 
running particle is not retrapped until Ib < Iq, and the current- voltage characteristic is 
hysteretic. In the presence of noise, the SQUID can be thermally excited out of the zero- 
voltage state at bias current Ib < Io [36] - In addition, the weak-link can macroscopically 
quantum tunnel through the barrier from the zero- voltage to the voltage state [37] . 

A Superconducting Quantum Interference Device (SQUID) is a superconducting 
ring interrupted by at least one Josephson weak link. For SQUID microscope 
applications a SQUID with two weak-links is most often used. Such a SQUID is labelled 
a dc-SQUID, since the modulation of the SQUID critical current with applied magnetic 
field can be observed using time-independent bias currents. Consider a SQUID with two 
weak-links 1,2 with critical currents I c i,I C 2, phase drops across the weak-links ip±,ip2, 
and inductances in the two arms of the SQUID L 1; L 2 . Neglecting for the moment the 
current through the resistive and capacitive elements of the junctions, the bias current 
across the SQUID Ib and the magnetic flux through the SQUID $ are given by 

Ib = h sin <pi + I 2 sin <p 2 

$ = $ a - L x h + L 2 J 2 , (4) 

where $ a is the externally applied flux. Integrating all the changes in phase around the 
SQUID loop along a path sufficiently deep inside the superconductor that there are no 
supercurrents along the integration path (i.e. they are shielded out by the bulk of the 
superconductor), and using the fact that the canonical momentum of the Cooper pairs 

is given by ^canonical = V + where q = 2e is the charge on the Cooper pairs and A 
is the vector potential, and finally by insisting that the Cooper pair wavefunction be 
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Figure 1. Plot of the critical current I c of a symmetric (Li = L2 = L/2, I i=i2=Io) 
2-weak link SQUID as a function of applied flux, for different values of the parameter 
7 = 27ri/o/ < £ > o- lo is the critical current of one of the weak-links, and L is the total 
inductance of the SQUID. 



single valued, leads to 

2 7m = ^ 2 -^ 1 + _($ a + L 2 i 2 -£ 1 /i) (5) 

n an integer. Solutions to Equations H] and are shown in Figure [TJ for the case of 
symmetric weak-link critical currents (I c \ = I&) and arm inductances (L\ = L2). 
The parameter 7 = 2nLI /^ (L^Lx+L^) determines the depth of the modulation 
of the SQUID critical current with flux: larger SQUID inductances lead to smaller 
modulations. The solutions for non-symmetric SQUIDs are more complicated, but 
qualitatively similar [38] to those shown here. 

There are 4 relevant sources of noise in dc SQUIDs: Johnson noise arises from 
thermal fluctuations, assumed to be dominated by the shunt resistor in the RSJ model 
[39J; shot noise arises from the discrete charge of the quasiparticles [40] and Cooper 
pairs [UJ traversing the weak-links; quantum noise arises from zero point motion in 
the shunt resistors [42J ; and 1/f noise can arise from a number of sources. Tesche and 
Clarke calculate [39] that the signal to noise ratio in dc SQUIDs is optimized when 
Pl = 2L/ /$o ~ 1- in addition, the Stewart-McCumber parameter (3 C should be about 
1, so that the SQUID is just non-hysteretic. If these two conditions are met, the optimal 
flux noise power S<s> from the first three sources is given by: 

f AkBTLinLC) 1 / 2 Johnson noise 
S§ = < hL Shot noise (6) 

I HL Quantum noise. 

Typical values for scanning SQUID sensors using a Nb-A^Os-Nb trilayer process are L 
= 100 pH, I = 6 11k, C = 8.36 x 10~ 13 F, R = 5 Q. and T = 4.2 K. These parameters 
lead to ^Johnson = v^Sj = 3 x 10 _7 $o/V / Hz for the flux noise due to Johnson noise; 
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$shot = 1-2 x 10 _7 $ /vHz for shot noise; and $ quantum = 5 x 10~ 8 $o/v#z for noise 
due to zero point motion. 

In the frequency range of interest Johnson noise, shot noise, and quantum noise 
are "white": the frequency distribution of the noise is independent of frequency /. 
1/f noise, as its name implies, has a noise amplitude that follows a 1/f dependence. 
It almost always dominates the noise in SQUIDs at sufficiently low frequencies. In 
the Dutta, Dimon, and Horn model |43j 1/f noise results from the superposition of 
"shot-like" contributions from two-level traps with a large distribution of trap energies. 
These traps could be, for example, defects in the oxide in a tunnel junction, or trap 
sites for superconducting vortices in the superconducting bulk making up the SQUID. 
Koch, Divincenzo, and Clarke jH] have suggested that 1/f noise in SQUIDs could result 
from coupling of flux into the SQUID from electronic spins with a wide distribution of 
characteristic times. 

As can be seen in Figure HJ the critical current of the SQUID is periodic in flux, with 
a period given by $o- The most commonly used scheme for determining this magnetic 
flux uses ac flux modulation in a flux locked loop jUJ HB]- If the SQUID is current 
biased just above the critical current the voltage varies sinusoidally with flux. If the dc 
component of the total (external+feedback) flux is such that the SQUID voltage is at an 
extremum, the first harmonic voltage response will be zero and any flux change will give 
a linear response. Negative feedback on the dc component of the flux is then used to keep 
the response zeroed: this feedback signal is linearly proportional to the flux threading 
the SQUID, with a constant of proportionality that can easily be determined with great 
precision. Traditionally the voltage signal from the SQUID has been amplified at low 
temperatures using a tuned LC circuit or a transformer, and at room temperature using 
phase sensitive lock-in detection. However, SQUID array amplifiers can also be used 
[17] HH] . These have the advantage of bandwidths of over 100 MHz, and no requirement 
for ac flux modulation. 

For this review I will concentrate on high-spatial resolution magnetic microscopies. 
High spatial resolution is achieved by placing a small sensor close to the sample. In the 
case of SQUID microscopy this means making the SQUID small, or making a pickup loop 
integrated into the SQUID small. The scaling of the flux signal and spatial resolution 
with sensor size depends on the type of field source. For this review I will concentrate 
on three basic field sources: A point dipole in generates the magnetic induction 

Ho 3f(f • m) — m 
Atc r 3 

where f is the direction and r is the magnitude of the vector between the dipole source 
and the point of interest. The magnetic flux $ s through a square area of side d oriented 
parallel to the xy plane and centered at r = zz above a dipole at the origin with moment 
m oriented parallel to the z axis is given by 

$ 5 = , 2/W * . Dipoleflux (8) 

7rJd 2 /2 + z 2 (d 2 + Az 2 ) 



B = ^ , Dipolefield (7) 



In the limit z — > $ s — > 2y/2fi m/rrd. 

A superconducting vortex generates the induction 

B = -r, Monopole field (9) 

Z7TH 

where $o = h/2e is the total magnetic flux generated by the vortex. The magnetic flux 
$ s through a square area of side d oriented parallel to the xy plane and centered at 
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Figure 2. Plot of the flux through a square area of side d a height z above a (a) dipolc 
source with moment m oriented normal to the scan plane (b) a monopole source with 
flux 4>o, and (c) a line source with current /. For a given value of z/d the flux signal 
scales like d _1 for a dipole source, is independent of d for a monopole source, and scales 
like d for a current line source. 



r = zz above a vortex at the origin is given by 

2$ ( d 2 \ 

$ s = — -tan" 1 . Monopole flux (10) 

7T \2zV2d 2 + Az 2 J K ' 

In the limit z — » $ s — > $o- Finally, an infinitely long and narrow line of current / in 
the y direction at x = 0, z = generates the field 

B = —9, Current field (11) 

2-KT ' V ' 

where 9 = (xz — zx)/r. The magnetic flux $ s through a square area of side d oriented 
parallel to the xy plane at a height z above this line has peaks at x = ±\/d 2 + z 2 /2 
with magnitude 



$ s = ^ln 



(d - Vd 2 + z 2 ) 2 + Az 2 



Current flux (12) 



Xd + Vd 2 + z 2 ) 2 + 4z 2 _ 

In the limit z — > the magnetic fluxes at the peaks diverge logarithmically. 

Figure [2] plots the calculated flux through a square area of side d, as a function of 
x, the lateral position of the SQUID relative to the field source, for various spacings z 
between pickup loop and sample, for these three different sources of field. For a given 
value of z/d for a dipole source the flux signal through the SQUID gets larger as d gets 
smaller, scaling like 1/d. For a monopole source, such as a superconducting vortex, the 
peak flux is independent of d at constant z/d. For a line of current, the peak flux signal 
increases like d as the pickup loop diameter increases. In all cases the spatial resolution 
is about d for z « d. 

Figure plots the ratio of the peak flux at a sensor height z to the peak flux at 
z/d = 0.1 for the three different field sources, while Figure Eb plots the full width at half- 
maximum Ax for the dipole and monopole field sources, and the peak-to-peak distance 
for the current line source, as a function of sensor height z. The peak amplitude falls 
off nearly exponentially with sensor height for a dipole source, but falls off less quickly 
for a monopole or current line source. Similarly, the spatial resolution falls off quickly 
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Figure 3. (a) Plot of the ratio of the flux through a square area of side d at a height 
z to that at a height z/d — 0.1 for various sources of field, (b) Full width at half- 
maximum for the dipole and monopole field sources, and the peak-to-peak distance 
for the current line source, of the flux through a square area of side d as a function of 
height z. 



with sensor height for the dipole source, but less quickly for the monopole and current 
line sources. 

The same qualitative conclusions can be made for more complex sources of field 
[49J. Consider a geometry in which the sample takes up a half-space z < 0. For the 
static case the magnetic fields in the half space z > can be written as the gradient 
of a scalar potential (p m : B = \/f m , where \/ 2 f m = 0. If we resolve the fields at the 
surface of the sample into Fourier components 

hk{z = 0) = (2^p J oo i-oo jB(x ' y ' z = °y {hxX+kyy)dx d v> ( 13 ) 

each of the components of the field for z > will be of the form bk{z) = bk(z = 0)e~ kz , 
where k = ^Jk% + ky\ The fields will decay exponentially as z gets larger, with the 
higher Fourier components decaying more rapidly. Therefore to get optimal sensitivity 
and spatial resolution, the sensor must be placed close to the sample. 

It was recognized almost immediately after the demonstration of the Josephson 
effect [50] and superconducting quantum interference effects [51] that local magnetic 
fields could be imaged by scanning samples and SQUIDs relative to each other 
[52J. However, it was not until 1983 that the first two-dimensional scanning SQUID 
microscope was built by Rogers and Bermon at IBM Research [53]. This instrument was 
used to image superconducting vortices in association with the IBM Josephson computer 
program. Other notable early efforts in SQUID microscopy were by the Wellstood group 
at the University of Maryland (5U [551 [56] , the van Harlingen group at the University of 
Illinois [571 EH], the Wikswo group at Vanderbilt University [59], the Clarke group at U.C. 
Berkeley [60] and the Kirtley group at IBM Research [611 [62]. Kirtley and Wikswo [16] 
reviewed some of fundamentals and early work concerning scanning SQUID microscopy. 

There are several competing strategies for achieving good spatial resolution in a 
SQUID microscope sensor. The first strategy is to make the SQUID very small (//- 
SQUID), with narrow and thin constrictions for the Josephson weak links [631 El]- 
This strategy has the advantage of simplicity, since only one level of lithography is 
required. The original yU-SQUIDs had hysteretic current-voltage characteristics. This 
meant that flux-locked-loop feedback schemes could not be used, and the flux sensitivity 



9 





Figure 4. Advanced integrated scanning SQUID susceptometer design, (a) Diagram 
of a counterwound SQUID susceptometer. (b) Optical micrograph of pickup loop 
region of an optically denned SQUID susceptometer. The outer ring is the field coil 
and the inner ring is the flux pickup loop of the SQUID. The edge of the silicon 
substrate is just visible on the left side. The inset shows an atomic force microscopy 
cross-section of the structure along a horizontal line through the center of the pickup 
loop. The pickup loop is closest to the sample when the tip is aligned at a 2.5° angle, 
(c) Design for a SQUID sensor in which the pickup loop is defined using focussed ion 
beam (FIB) etching. In this case the pickup loop can touch down first if the alignment 
angle is between 2° and 5°. (d) Scanning electron microscopy images of the sensor 
after FIB definition of the pickup loop. Reprinted figure with permission from N.C. 
Koshnick, M.E. Huber, J. A. Bert, C.W . Hicks, J. Large, H. Edwards and K.A. Moler, 
|Appl Phys. Lett. 93, 243101 (2008)] Copyright 2008 by the American Institute of 
Physics. 



of these SQUIDs was relatively poor. However, recently non-hysteretic, sub-micron sized 
SQUIDs have been made [65j ESI E7J . Hao et. al [67] made their SQUIDs non-hysteretic 
by using a second layer of tungsten to shunt them, and reported white noise floor levels 
of O^/^o/Hz 1 / 2 in a SQUID with a diameter of about 370 nm. 

A second strategy is a self-aligned SQUID recently reported by Finkler et al. [68J. 
In this work three aluminum evaporations are made onto a quartz tube that has been 
pulled into a sharp tip with apex diameter between 100 and 400 nm. The first two 
evaporations, performed at an angle of 100 degrees relative to the axis of the tube, form 
the leads. A third evaporation along the tip axis forms a ring with two weak links 
to the leads, forming a dc SQUID. Finkler et al. report non-hysteretic current- volt age 
characteristics and a flux sensitivity of 1.8x lO^c&o/Hz 1 / 2 for SQUIDs with an effective 
area of 0.34 /zm 2 , operating at fields up to 0.6 T. 

A third strategy is to make a more conventional SQUID, but to have well shielded 
superconducting leads to a small pickup loop integrated into it [62JEnjE2] ( see Fig- BJ- 
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This has the disadvantage of complexity, since multiple levels of metal are required to 
shield the pickup loop leads, but has the advantages of reduced interaction between the 
SQUID and the sample and ease of using flux feedback schemes. Figure H] illustrates 
the properties of current advanced scanning SQUID sensor design using the integrated 
pickup loop strategy [71] . This SQUID sensor uses a Nb-A10 x -Nb trilayer for the 
junctions, two levels of Nb for wiring and shielding, and SiC>2 for insulation between the 
various levels, as well as local field coils integrated into the pickup loop region, to allow 
local magnetic susceptibility measurements [69l 172] . High symmetry in the SQUID and 
field coil layout (see Fig. UK) and tapered terminations are desirable to discriminate 
against background fields and reduce unwanted electromagnetic resonances [75]. Fig. 
Hb is an optical micrograph of a scanning SQUID susceptometer in which both the field 
coil and pickup loop were defined using optical lithography. Fig. Hb shows the layout for 
a device in which the pickup loop area is left as a "tab" in the optical lithography step, 
and then patterned into a loop using focussed ion beam lithography |74j . The completed 
device is shown in a scanning electron micrograph in Fig. HH. This device uses a non- 
planarized process, meaning that each successive layer conforms to the topography of 
the previous one. This leads to constraints on the alignment angles that can be used to 
get the minimum spacing between pickup loop and sample, as illustrated in the insets in 
Fig. Hb,c. These constraints can be reduced by using a planarized process [75], in which 
the insulating layers are chemically-mechanically polished to planes before successive 
Nb wiring levels are added. 

If we assume a SQUID or SQUID pickup loop diameter of l/mi, a SQUID-sample 
spacing of 0.1/im, and a SQUID noise level of 10 _6 $o/Hz 1 / 2 , and using the peak signals 
from Fig. [21 we find a minimum detectable dipole signal of 139 electron spins/Hz 1//2 , a 
minimum flux signal of 1.25xlO- 6 <£> / Hzl/2 ; and a minimum detectable line of current 
signal of 1.65 nA/Hz 1 / 2 . 

2.2. Hall bar microscopy 

A Hall bar develops a transverse voltage V x = —I y B z /n2d£ when a current I y passes 
through it in the presence of a magnetic induction B z perpendicular to the plane 
of the Hall bar [30| |3T| 132] , where ri2d is the carrier density per unit area of the 
Hall bar. Therefore materials with small carrier densities such as semi-metals, semi- 
conductors, or two-dimensional electron gases at the interface between semi-conductors 
with different bandgaps develop larger Hall voltages. Early scanning Hall bar systems 
used evaporated films of bismuth [76], InSb [77J, or GaAs [78]. More recently 
high sensitivity and spatial resolution have been achieved with GaAs/Al^Gai-^As 
hetero junction structures JT9J EQl [8TJ [82j |83j El]. In addition, 250 nm sized scanning 
probes using GaSb/InAs/GaSb [83], and micron-scale Si/SiGe and InGaAs/InP Hall 
crosses have been characterized at low temperatures [SSI E5]- Hicks et al. [S7J have 
fabricated Hall bars as small as 85 nm on a side (see Fig. [5]) from GaAs/ALGai-aAs 
hetero junction electron gas material, and estimate a field noise of 500/iT/Hz 1 / 2 and a 
spin sensitivity of 1.2 xIO^hb/ Hz 1 / 2 at 3 Hz and 9K for sensors 100 nm on a side. Sandhu 
et al. have fabricated 50 nm bismuth Hall bar sensors with a noise of 0.8G / Hz 1 ! 2 |88j. 
At present Hall bars from GaAs/Al^Gai-^As hetero junctions are less sensitive above 
about 100 K because of thermal excitations. Scanning Hall bar sensors made of InAs 
have been used to image single magnetic biomolecular labels at room temperature [89J . 

Magnetic tunnel junctions and Hall bars in the ballistic regime [901 |9TJ |92] produce 
a signal which is proportional to < B z >, the magnetic field perpendicular to the 
sensor, averaged over the sensor. The situation is more complicated for a Hall bar in 
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Figure 5. Scanning electron microscope image of a 100 nm Hall probe. The four leads 
are separated by narrow etch lines. The gate shields the Hall cross from stray electrical 
charges and allows the modulation of the 2DEG beneath. The Hall probe will touch the 
sample surface at the contact tip. Reprinted figure with permission from C.W. Hicks, L. 



Luan, K.A. Moler, E. Zeldov, and H. Shtrikman, Appl. Phys. Lett. 90, 133512 (2007) 



Copyright 2007 by the American Institute of Physics. 



the diffusive regime [93], but may be qualitatively similar. The average fields < B z > for 
the three sources of field can be inferred from Eq.slSl fl0land[T2l (Fig. [2]), by substituting 
$ s — >< B z > d 2 and in Fig. [3] by substituting $ p — >< B z > p , where < B z > p is the 
peak value of the field averaged over the sensor area. For a given value of z/d, < B z > 
is proportional to 1/d 3 for a dipole source, proportional to 1/d 2 for a monopole source, 
and proportional to 1/d for a line of current source. Boero et al. [91] estimate that the 
optimal Johnson noise limited minimum detectable magnetic field in a micro-Hall cross 
geometry is given by 

B mm ~ Thermal (14) 

VsatW 

where fee is Boltzman's constant, T is the operating temperature, Rq is the output 
resistance at zero magnetic field, v sa t is the saturation carrier drift velocity, and w is 
the width of the cross. This results in B min ~ 2nT / y/Wz for w=l/im at 300K for doped 
InSb. Thermal noise is expected to be independent of sensor size, and therefore the 
scaling for sensitivity is the same as for < B z >. At the low temperatures used for 
imaging superconductors, transport in the Hall sensor is often quasi-ballistic and the 
signal-to- noise is better characterized by a mobility and an optimum probe current. The 
latter is typically defined as the maximum current before the onset of pronounced 1/f 
noise due to one of several possible sources, e.g., generation- recombination noise at deep 
traps or across heterointerfaces, heating or impact ionization. However, noise in micron 
scale Hall bars is often dominated by 1 // noise EZ] • The minimum detectable field 
can then be written as 189 



1 a H G N Af 
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Figure 6. MFM modeling using the point dipole tip approximation: Force derivative 
curves for (a) a dipole source located at the origin with magnetic moment m s oriented 
parallel to z (solid line) or parallel to x (dashed line), Eq. (b) a monopole source 
located at the origin with magnetic flux $o, Eq. [9j and (c) a line of current / along 
the y axis, Eq. 111! The tip is modeled as a point dipole with moment m* oriented 
parallel to the z-axis scanning a height z above the xy plane along the x axis (Eqn.'s 

OH). 



where /x is the mobility, d is the sensor size, an is Hooge's 1/f noise parameter, 
Gjy ~ 0.325 is a constant, n is the carrier density, Af is the measurement bandwidth, 
and / is the frequency. Eq. [15] combined with Eq.'slH [101 and [T2l imply that for a Hall 
bar dominated by 1/f noise, for a given z/d the minimum detectable dipole moment 
scales like d 2 , the minimum detectable flux scales like d, and the minimum detectable 
current is independent of d. It has been reported that low frequency noise in Hall bar 
devices can be significantly reduced by optimizing the voltage on a gate over the Hall 
cross [22 EZ] and that sufficiently small devices have noise which is composed of a single 
Lorentzian spectrum [96] . 

2. 3. Magnetic force microscopy 

A magnetic force microscope [971 EH [99] images small changes in the resonance frequency 
of a cantilever due to the interaction between magnetic material on a sharp tip at the 
end of the cantilever and local sample magnetic fields. The signal from a magnetic force 
microscope is proportional to dF tiPtZ /dz, the derivative of the z-component of the force 
on the tip with respect to its z-position. Since this force sums contributions from all the 
magnetic material in the tip, the sensitivity and spatial resolution of MFM depends on 
the shape of the tip. For some applications and tip materials the tip can be modeled as 
having a magnetic monopole at its apex, with the force approximated by F tiPiZ = qB z , 
where q is the effective magnetic monopole moment of the tip [100[ [32] . However, Lohau 
et al. have shown that their experimental MFM data can be modeled either using a 
monopole or a dipole approximation for the end of the tip, but that in either case the 
effective position of the monopole or dipole is a function of the field gradient from the 
sample jlOlj . 

It is instructive to model the ultimate spatial resolution and sensitivity of an MFM 
in the point monopole and dipole approximations for the tip. The z-component of the 
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force derivative on the tip from a field source can be written as [101J 

+ m t,x + m^y^T^r + m^z^T^, (16) 



dF z _ q dB z d 2 B x d 2 B y d 2 B z 

—7, — — ; 1" TTlt x _ „ — h m t y _ 5 — h m t z 

az /i az as 2 az^ az z 



where g is the magnetic monopole flux, and mt t i are the dipole moments of the tip in 
the x, y, z directions. Assume for simplicity that the tip is magnetized only in the z- 
direction. Then m t , x = m t y = and the force derivatives due to the dipole, monopole, 
and current line sources of Eq. 's [7J [91 and [11] are given by 

dF z 3m tz Ho 2 2 

~dz = 4nr 9 ^ ~ 5 ( m ^' x + m S:y y)z{3r - 7z ) 

+ m SjZ (3r 4 — 30r 2 z 2 + 35z 4 ) Dipole — dipole 

3$o m t z z 5z 2 — 3r 2 



2tt r 7 
m^z/iolx x 2 — 3z 2 



Dipole — Monopole 
Dipole — Current line (17) 



7T [x 2 + Z 2 ) 3 

where m tj2 is the dipole moment of the tip. 

Similarly, the expressions for the force derivative in the monopole approximation 
are given by: 

9F Z 3^ r , \t 2 r 2\ 

~dz~ = ~ W [ ^ Xms ' x + yms ^ r ~ 5z ) 

+ m SiZ z(3r 2 — 5z 2 ) Monopole — dipole 

2 oJl 



g$o r 2 - 3z 
2nfi r 5 
ql xz 



Monopole — Monopole 
, ; Monopole — Current line (18) 

71 [X 2 + Z 2 ) 2 

The force derivative relations given by Eq.'s[T7]and[T8]are plotted in Figures [6] and 
UJ respectively. 

Within these approximations the force derivative signal strength increases strongly 
with decreasing z - like 1/z 5 for a dipole source, l/z 4 for a monopole source, and 1/z 3 
for a current line source, in the point dipole tip approximation. In the monopole tip 
approximation the powers are 1/z 4 , 1/z 3 , and 1/z 2 for the three sources, respectively. 
The widths of the predicted force derivative features are proportional to the height z 
of the tip above the source. The constant of proportionality does not depend strongly 
on source: the full width at half-maximum (or distance between the maximum and 
minimum force derivative for a line of current source) is Ax/ z=0.6, 0.66, and 0.74 
for the dipole tip, and 0.74, 1, and 1.16 for the monopole tip, for the three sources 
respectively. Of course, these conclusions would be modified for a real tip with a finite 
size. It has been estimated that the effective tip-sample spacing and radius of curvature 
for conventional MFM tips are both about 10 nm [100] . 

In order to estimate the sensitivity of MFM, one needs to know the minimum 
detectable force derivative, which depends on a number of characteristics of the MFM, 
including those of the cantilever. Using the lumped mass model for an MFM cantilever, 
the equation of motion of the tip end position z is that of a driven, damped harmonic 
oscillator 

cPz dz 

«i— + T— + kz = F signa i(t) + F therma i(t), (19) 
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Figure 7. MFM modeling using the monopole tip approximation: MFM force 
derivative curves for (a) a dipole source located at the origin with magnetic moment m s 
oriented parallel to z (solid line) or parallel to x (dashed line), Eq. (b) a monopole 
source located at the origin with magnetic flux $o, Eq. [9J and (c) a line of current / 
along the y axis, Eq. 111! The tip is modeled as a point monopole with magnetic flux 
q scanning along the x axis a height z above the xy plane (Earns [T8|) . 



where m is the effective mass of the cantilever, V = k/u^Q is the damping constant, k 
is the cantilever spring constant, co is the angular resonant frequency, Q is the quality 
factor, F signal is the signal force (in this case the magnetic interaction between the tip 
and sample), and -Fthermai is the force due to thermal fluctuations. The equipartition 
theorem implies that the power spectral density of thermal fluctuations Sf = ^TkuT. 
Then the thermally limited minimum detectable force derivative is given by 

dF = 1 j 4kk B TBW 

dz min 

where BW is the measurement bandwidth and A is the amplitude of oscillation of the 
cantilever. 

A major thrust in research in MFM is to reduce the spatial extent of the magnetic 
material in the cantilever tip to improve spatial resolution and spin sensitivity. This is 
often done by shaping previously deposited material using focussed ion beam etching 
[102[ I103j . Sharply defined magnetic tips can also be produced by evaporating magnetic 
material on an electron beam deposited carbon and ion etched needle |1U4] . or onto a 
carbon nanotube |105] at the end of a conventional Si cantilever tip (see Figure[8]). Small 
amounts of magnetic material can also be deposited directly onto cantilever tips through 
nanoscale holes fabricated in a stencil mask |1U6] . One approach would be to place nano- 
magnets directly on the apex of cantilever tips. This would serve three goals: improved 
spatial resolution, improved spin sensitivity, and improved interpretability, since the 
magnetic moment of the deposited nano-magnets would have a narrow distribution 
[107] . However, such tips could be especially susceptible to switching of the tip moment 
in sufficiently strong magnetic fields, reversing the contrast of the MFM images [108] . 
Although this complicates the interpretation of the images, the tip magnetization state 
is also known better, since the tip is generally close to saturation either up or down 
except for a narrow window around the switching points. The tip in MFM exerts 
a relatively strong magnetic field, comparable to the saturation magnetization of the 
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Figure 8. Scanning electron microscope (SEM) image of a MFM cantilever. Insert 
upper left: high-resolution SEM image of the apex of the pyramid, where a coated 
carbon nanotube tip (CCNT) is visible. The arrow shows the direction of the metal 
evaporation. Reprinted figure with permission from Z. Deng, E. Yenilmez, J. Leu, J.E. 
Hoffman, E.W.J. Straver, H. Dai, and K.A. Moler, |Appl. Phys. Lett. 85, 6263 (2004) 
Copyright 2004 by the American Institute of Physics. 

magnetic material used in the tip, on the sample. This can cause switching of the local 
magnetic moment of the sample |109j . but can also be taken advantage of for vortex 
manipulation studies (see Section [375]) . 

Assume for the moment that it will be possible to locate small nanoparticles at the 
end of an MFM cantilever tip. A 7 nm diameter cobalt nanoparticle will have a total 
dipole moment of m t = 2.5 x 10~ 19 v4 — m 2 , assuming a magnetization of M — 1.4 x 10 6 
A/m pIO]. Taking values of C?=50,000, w =27rx50 kHz, A=l nm, A;=2N/m, and 
T=4K, the thermally limited force gradient is dF/dz\ m i n =1.7x 10 _7 N/m. Using the 
dipole approximation for the tip, the peak force due to a sample dipole m s with both 
tip and sample moments oriented parallel to the z axis is given by (Eq. [T7J Fig. [6]) 
dF/dzp ea k = 24fi m s m t /4:7rz 5 . Using z=15 nm results in a minimum detectible sample 
moment of 23 Bohr magnetons. From Fig. [6] the full width at half maximum of the 
force derivative response would be Ax = 0.6z = 9nm. 

2.4- Summary 

Each of the three techniques described here has advantages and disadvantages: MFM 
has the best spatial resolution (~10-100 nm), SSM the worst (~ 0.3-10/im), with Hall 
bars intermediate (~ 0.1-5/xm). The relative sensitivities of the three techniques depend 
on the type of field source. As in all microscopies there are tradeoffs between sensitivity 
and spatial resolution. The spatial resolution for SQUIDs and Hall bars is roughly given 
by Ax ~ d for heights z « d, where d is the sensor size (Fig.'sO [3]). For MFM the 
spatial resolution is roughly set by the tip-sample spacing z, in the limit where the 
radius of curvature of the tip end is smaller than z (Fig.'s [7f6|) . The flux sensitivity 
of SQUIDs is nearly independent of SQUID or pickup loop size, so the field sensitivity 
is roughly proportional to the inverse of the pickup loop or SQUID area. With Hall 
bars limited by thermal noise the field sensitivity is independent of sensor size, but Hall 
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bars dominated by 1/f noise are expected to have a field sensitivity that scales like 1/d. 
Current MFM's measured field gradients, rather than fields, and the tradeoff between 
sensitivity and spatial resolution is less clear than for SSM and SHM, although larger 
amounts of magnetic material will produce larger sensitivity but will require a larger 
tip volume, reducing resolution. Table [1] compares sensitivity and spatial resolution for 
a state of the art SQUID, with a pickup loop 0.6/im in diameter and a flux noise of 
0.7xlO- 6 $o/-^<z 1,/2 (Figure H] [TT] ) , a state of the art Hall bar with sensor size 100 nm 
and a field noise of 500/iT/ifz 1/2 at 3 Hz (Fig. E [87J, and a hypothetical MFM with a 
7 nm Co nanoparticle at the tip as discussed in Section 12.31 These numbers should be 
treated with a great deal of caution. For example, I have assumed that z/d = 0.1 for the 
SQUID and SHM, and z=15 nm for the MFM. These may be unreasonably optimistic 
estimates. However, the Table [T] does indicate how relative sensitivities depend on the 
field source. As an example, SSM and MFM have roughly comparable spin sensitivities, 
but SSM has nearly two orders of magnitude better sensitivity for a monopole source, 
and nearly 4 orders of magnitude better sensitivity for a line of current source. 



Table 1. Sensitivities and spatial resolutions 





Dipole 




Monopole 


Line of current 




N(/jl b / VHz) 


Ax(nm) 


$(<J>o/\/#I) Ai(nm) 


l(A/VJTz) Ax(nm) 


MFM 


23 


9 


3.4xl0" 5 10 


2.6xl0" 5 11 


SHM 


4.5xl0 4 


100 


3xl0- 3 108 


1.14xl0~ 4 118 


SSM 


78 


600 


8.8xl0~ 7 650 


5.5xl0~ 9 710 



Table [2] provides a summary of the scaling exponents n for the minimum detectable 
field sources for SSM, SHM, and MFM. In this table the minimum detectable dipole 
moment, flux, and current are proportional to d n for SSM and SHM, while for MFM 
these quantities are proportional to z n . For SSM and SHM it is assumed that z/d is 
held constant. For SSM it is assumed that thermal noise dominates. There are two 
columns for SHM, for noise dominated by thermal fluctuations, or by 1/f noise. For 
MFM there are two columns, for the monopole and the dipole tip approximations. 



Table 2. Scaling exponents n 
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There are further considerations in deciding which magnetic microscopy is best for 
a given application: SSM requires a cooled sensor, about 9K for the most technically 
advanced Nb SQUIDs, although Hall bars and MFM also have the highest sensitivities 
at low temperatures. SSM is the most straightforward to calibrate absolutely. Hall bars 
tend to be mechanically delicate and the most sensitive to destruction by electrostatic 
discharge. MFM applies the largest fields to the sample. In general Hall bars and MFM 
can tolerate larger applied magnetic fields than SQUIDs. As we shall see in Section [3j 
each has been very useful in the study of superconductors. 
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Figure 9. (a) Magnetic force microscope image of 0.3 /mi x 0.3 jum areas for a Nb 
film and the corresponding topographic image (b). (d) MFM image of a single vortex 
(2x2 /j,m 2 ) in a YBCO film and (e) the corresponding topographic image. The white 
circles in (b) and (e) correspond to the distance at which the stray field emanating 
from the vortex has decreased to 1/e of its maximum value. The small dotted circle 
in (b) defines an area with diameter 2£ ps 22 nm. The lower panels (c) and (f) are 
cross-sections of the MFM data and the topography profiles along the white lines 
indicated in (a), (b), (d), and (e). The dotted curves correspond to a Gaussian profile 
fitting of the MFM signal after filtering with a low-pass filter. Reprinted figure from 
Physica C: Superconductivity 369, A. Volodin, K. Temst, Y. Bruynseraede, C. van 
Haesendonck, M.I. Montero, I.K. Schuller, B. Dam, J.M. Huijbregtse and R. Griessen, 
"Magnetic force microscopy of vortex pinning at grain boundaries in superconducting 
thin films", P. 165-170 (2002), with permission from Elsevier. 



3. Applications 

3.1. Static imaging of magnetic flux 

Understanding the trapping of vortices in superconductors is of crucial technological 
importance, both because vortex pinning is the primary mechanism for enhanced 
critical currents in type II superconductors [111] , and also because trapped vortices are 
an important source of noise in superconducting electronic devices. Superconducting 
vortices are especially easy to image using magnetic imaging because they are highly 
localized, quantized, and have relatively large magnetic fields. A number of different 
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Figure 10. SQUID microscope images of 35 /im wide YBCO thin film strips cooled 
in magnetic inductions of (a) 5, (b) 10, (c) 20, and (d) 50 /u,T. e) Critical inductions 
for vortex trapping in YBCO thin films as a function of strip width. The squares 
(B c _|_) represent the lowest induction in which vortices were observed, and the dots 
(B c _) are the highest inductions for which vortex trapping was not observed. The 
dashed-dotted line is the metastable critical induction Bo [Eq. [23] (1 20] . the short- 
dashed and long-dashed lines are B^ [Eq. [22], the critical induction calculated using 
an absolute stability criterion, and the solid line is Bk [Eq. [24], calculated using a 
dynamic equilibrium criterion between thermally activation and escape for vortices. 
Reprinted figure with permission from K.H. Kuit, J.R. Kirtley, W. van der Veur, C.G. 
Molenaar, F.J.G. Roesthuis, A. CP. Troeman, J.R. Clem, H. Hilgenkamp, H. Rogalla, 
and J. Flokstra, |Phys. Rev. B 77, 134504 (2008)[ Copyright 2008 by the American 
Physical Society. 



techniques have been used for imaging superconducting vortices, including SSM, SHM, 
and MFM. MFM has the advantage for this application that it can image both the vortex 
magnetic fields and sample topography simultaneously, so that macroscopic pinning 
sites can be identified. It has been used for a number of years to image superconducting 
vortices pZE2l E31 E3 ESI ESI EH E3 El]. Care must be taken however, because 
the large magnetic fields exerted by the tip can distort and dislodge vortices. As we will 
see in Section 13.51 this can also be an advantage. An particularly striking example of 
simultaneous magnetic and topographic imaging of vortices in Nb and YBCO is shown 
in Figure [9] [EE]. 

3.1.1. Narrow strips The sensitivity of high-T c superconducting sensors such as 
SQUIDs and hybrid magnetometers based on high-T c flux concentrators is limited by 
1/f noise. One source of this noise is the movement of vortices trapped in the sensor. 
One technique for reducing this source of noise is to trap vortices at a distance from 
the sensitive regions of the superconducting circuitry using holes or motes [121[ 1122] . 
Jeffery et al. [123j used a scanning SQUID microscope to image flux trapping in 
superconducting electronic circuitry with motes. Veauvy et al. |124j studied vortex 
trapping in regular arrays of nanoholes in superconducting aluminum. Noise from 
trapped vortices in high-T c SQUID washers and flux concentrators can be eliminated by 
dividing the high-T c body into thin strips |125[I126] . Below a critical magnetic induction 
no vortex trapping occurs in these strips for a given strip width. A number of models for 
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the critical inductions of thin-film strips have been proposed |127[ll20lll28[ll29] . Indirect 
experimental testing of these models was done by observing noise in high-T c SQUIDs 
as a function of strip width and induction |125[ I126[ 1130] . More direct experimental 
verification of these models was presented by Stan et al. [131] using scanning Hall 
probe microscopy on Nb strips and Suzuki et al. [132J using SQUID microscopy on 
NdBa2Cu30j / thin-film patterns with slots. Both experiment and theory found that 
the critical magnetic induction varied roughly as 1/W 2 , where W is the strip width. 
However, the experimental |131] and theoretical |12Q[ 1127] pre-factors multiplying this 
1/W 2 dependence differ significantly. 

Figure [TUa-d shows SQUID microscope images of 35/xm wide YBCO strips cooled in 
various magnetic inductions. Fig. [TUfe plots B c+ , the lowest induction at which vortices 
are trapped in the strips (squares), and B c _, the highest induction at which vortices are 
not trapped (dots) as a function of strip width. The lines represent various theoretical 
predictions. 

The Gibbs free energy for a vortex in a superconducting strip of width W in an 
applied induction B a perpendicular to the strip plane can be written as [133]: 

\aW . fnx\] <S> (B a -n$ 
&{x) = -In — —sin 



27T/ioA 



=F a ' *(W-x), (21) 



where A = 2X 2 /d is the Pearl length of the strip, A is the London penetration depth, 
a is a constant of order 1, n is the areal density of vortices, £ is the coherence length, 
and x is the lateral position of the vortex in the strip. The critical induction model 
of Likharev |127j states that in order to trap a vortex in a strip the vortex should be 
absolutely stable: the critical induction then happens when the Gibbs free energy in the 
middle of the strip is equal to zero, and leads to 

Bl = 2*Lu(!M). (22) 

This is plotted as the dashed and dotted lines in Fig. [10] for two assumed values of a. 
A second model, proposed by Clem |120] considers a metastable condition, in which the 
induction is just large enough to cause a minimum in the Gibbs energy at the center of 
the strip, d 2 G{W/2)/dx 2 = 0, leading to 

This is plotted as the dot-dashed line in Fig. [TUJ Finally, Kuit et al. proposed that 
the critical induction should result from a dynamic equilibrium between vortex thermal 
generation and escape. This leads to the critical induction 

Bk = 1-65^. (24) 

It can be seen in Figure \IU\ that the dynamic equilibrium model (solid line) fits the 
data the best of the three models, with no fitting parameters. Kuit et al. also found 
that n, the areal density of trapped vortices, increased linearly with applied induction 
above the critical value, in agreement with the dynamic equilibrium model, and that the 
trapping positions showed lateral ordering, with a critical field for formation of a second 
row that was consistent with a numerical prediction of Bronson et al. |134] . Similar 
good agreement between experiment and the dynamic equilibrium model was found by 
the same authors for Nb thin film strips. 
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Figure 11. Top row: Scanning Hall probe microscopy images of vortex configurations 
in square areas f 9.6 /im on a side in a superconducting Nb wire grid of 0.95 fim x 0.95 
/an square holes, at the filling fractions / indicated. The second row shows maps of the 
positions in the grid occupied by vortices. Note the ordering observed near / = 1/3 (c), 
/ = 1/2 (e), and / = 4/3 (h). Reprinted figure with permission from H.D. Hallen, R. 
Seshadri, A.M. Chang, R.E. Miller, L.N. Pfeiffer, K.W. West, C.A. Murray, and H.F. 
Hess, |Phys. Rev. Lett. 71, 3007 (1993)[ Copyright 1993 by the American Physical 
Society. 



3.1.2. Superconducting wire lattices, clusters, and nano-hole arrays Two dimensional 
arrays of Josephson junctions [1351 11361 11371 1138] 1139] fTlO] have been studied extensively, 
not only because they offer the possibility of studying Kosterlitz-Thouless effects |141] 
in an ordered 2- dimensional system, but also because of interesting effects that were 
predicted |142j and observed [1351 H36j to occur when the applied field is a rational 
fraction of the field required to populate each cell with a quantum of magnetic flux. 
Scanning magnetic microscopy has been used extensively for studying vortex trapping 
in two-dimensional superconducting wire networks [1431 158] . clusters [58], and ring 
arrays [1441 1145] . Superconducting wire networks are of interest because they are a 
physical realization of a frustrated xy model [1461 11421 1147] . where the variable is the 
phase of the superconducting order parameter. The filling fraction / of vortices/site 
can be varied by cooling the network in various fields, and detailed predictions have 
been made for the ground state vortex configurations at rational / fractions [148] . 
Studies of these systems were at first performed using transport [1491 11501 11511 1152] , but 
direct imaging of the vortex positions provides more detailed information. Runge et al. 
did magnetic decoration experiments on superconducting wire lattices [153j . Although 
these experiments provided direct information on the vortex positions, only one set 
of conditions, such as magnetic field and temperature, can be explored for each sample 
using this technique. An example of scanning Hall bar microscopy on a superconducting 
wire network is shown in Figure [TT1 [143] . In this figure ordered vortex arrangements can 
be seen near the fractional filling factors of /=l/3, 1/2, and 4/3, as predicted by Teitel 
and Jayaprakash [146] 1142] , along with grain boundaries between ordered domains and 
vacancies. 

Field et al. [154] studied vortex trapping in Nb films with large arrays of 0.3/mi 
holes on a square lattice. They found strong matching effects when the cooling fields 
corresponded to one or two vortices per site, as well as at several fractional multiples of 
the matching field. They also observed striking domain structure and grain boundaries 
between the domains. 

3.1.3. Half-integer flux quantum effect The discovery in 1986 of superconductivity 
at high temperatures in the cuprate perovskites |2] E] generated enormous excitement. 
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Figure 12. (a) Schematic diagram for the tricrystal (100) SrTiO*3 substrate used 
in the phase sensitive pairing symmetry tests of Tsuei et al. [61 j - Four epitaxial 
YBa2Cu307_,j thin film rings are interrupted by 0, 2, or 3 grain boundary Josephson 
weak links, (b) Scanning SQUID microscopy image of the four superconducting rings 
in (a), cooled in a field <5 mG. The central ring has $o/2 magnetic flux, where $o is 
the superconducting flux quantum, spontaneously generated in it. The other 3 rings, 
which have no spontaneously generated flux, are visible through a small change in 
the self-inductance of the SQUID when passing over the superconducting walls of the 
rings. 

Although there was a great deal of evidence that the superconducting gap in the cuprates 
was highly anisotropic, with a significant density of states at low energies |155l 1156] . 
phase sensitive measurements |157[ 1158] were required to distinguish, for example, d- 
wave from highly anisotropic s-wave pairing symmetry. The first such experiments 
|1591 1160] used the magnetic field dependence of the critical current in single junctions 
and 2-junction SQUIDs between single crystals of YE^CusOy-^ (YBCO) and Pb to 
demonstrate a 7r-phase shift between the component of the superconducting pairing 
order parameter perpendicular to adjacent in-plane crystal faces (e.g., ignoring the 
effects of twinning, between the a vs b axis normal faces) of the YBCO. Similar phase 
sensitive Josephson interference experiments were performed by Brawner et al. |161] , 
Iguchi et al. |162j . Miller et al. |163] and Kouznetsov et al. |164| . 

A phase-sensitive technique for determining pairing symmetry that is complemen- 
tary to Josephson interference is to image magnetic fields generated by spontaneous 
supercurrents in superconducting rings that have an odd number of intrinsic 7r-phase 
shifts in circling them, or in Josephson junctions that have intrinsic 7r-shifts along them 
[165[ H66] . The first such experiments were done by Tseui et al. [61], using YBCO 
films grown epitaxially on tricrystal SrTi03 substrates. A scanning SQUID microscope 
was used to image the magnetic fields generated by supercurrents circulating rings pho- 
tolithographically patterned in the YBCO films. In the original experiments the central 
ring (Fig. IT2l) . which has 3 grain boundary Josephson junctions, has either 1 or 3 intrin- 
sic 7r-phase shifts for a superconductor with predominantly d x 2_ y 2 pairing symmetry, 
and should therefore have states of local energy minima, when no external magnetic 
field is applied, with the total flux through the ring of $ = (n + l/2)$o, n an integer, 
for sufficiently large LI C products, where L is the inductance of the ring, and I c is the 
smallest critical current of the junctions in the ring. This is the half- integer flux quan- 
tum effect. The two outer rings that cross the tricrystal grain boundaries have or 2 
intrinsic 7r-phase shifts for a d x 2_ y 2 superconductor, and should therefore have local en- 
ergy minima at $ = n$>o - conventional flux quantization. Similarly, the ring that does 
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Figure 13. Scanning SQUID microscope image of an area of 1024 x 256 /im 2 along 
an asymmetric 45° [001] tilt YBa2Cu30y_5 bicrystal grain boundary (arrows). There 
are some ten bulk vortices in the grains (film thickness ^180 nm), but there is also 
flux of both signs spontaneously generated in the grain boundary. The bottom section 
shows a cross section through the data along the grain boundary measured in units 
of $o penetrating the SQUID pickup loop. Reprinted figure with permission from J. 
Mannhart, H. Hilgenkamp, B. Mayer, Ch. Gerger, J.R. Kirtley, K.A. Moler, and M. 
Sigrist, |Phys. Rev. Lett. 77, 2782 (1996)| Copyright 1996 by the American Physical 
Society. 



not cross any grain boundaries has no intrinsic phase drop, and also has local minimum 
energies for $ = n$o- These expectations were confirmed using a SQUID microscope 
using a number of techniques [STJ, consistent with YBCO having d x 2_ y 2 pairing sym- 
metry. Subsequent tricrystal experiments eliminated the possibility of a non-symmetry 
related origin for the half-integer flux quantization in the 3-junction ring |167j . ruled 
out extended s-wave symmetry |168] , demonstrated the existence of half-flux quantum 
Josephson vortices |169j . showed results consistent with d x 2_ y 2 pairing symmetry for 
the hole-doped high-T c cuprate superconductors Tl 2 Ba2Cu06+,5 |170j . Bi2Sr 2 CaCu208+<5 
|171j . and La 2 _ a ;Sr a; Cu04_j / |172] , for a range of doping concentrations |172] . and the 
electron-doped cuprates Ndi.gsCeo.isCuO^ (NCCO) and Prx .ssCeo.isCuO^ (PCCO) 
|173j . and for a large range in temperatures |174j . SQUID magnetometry experiments 
on two-junction YBCO-Pb thin film SQUIDs were performed by Mathai et al. |175j . and 
on tricrystals by Sugimoto et al. [176] . Josephson interference, SQUID magnetometry, 
and microwave measurements were made on biepitaxial all high-T c YBCO junctions by 
Cedergren et al. |177j . 

Grain boundary junctions in the high-T c cuprate superconductors are Josephson 
weak links. As such they are used for fundamental and device applications. In addition, 
grain boundaries play an important role in limiting the critical current density in 
superconducting cables made from the cuprates |178j . Anomalies in the dependence 
on magnetic field of the critical current of cuprate grain boundary junctions with 
misorientation angles near 45° can be understood in terms of their d x 2_ y 2 pairing 
symmetry combined with faceting, which naturally occurs on a 10-100nm scale in 
these grain boundaries |179j . This combination produces a series of intrinsic 7r-phase 
shifts along such grain boundaries. These 7r-shifts have spontaneous supercurrents and 
magnetic fields associated with them |18U] . An example of SQUID microscope imaging 
of these spontaneous fields is shown in Figure [13j If the distance between intrinsic tt- 
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Figure 14. Generation of half-flux quanta in connected and unconnected 
YBa2Cu307^.5-Au-Nb zigzag structures. Shown are scanning SQUID micrographs 
of (a) 16 antifcrromagnetically ordered half-integer flux quanta at the corners of a 
connected zigzag structure, and (b) 16 ferromagnetically ordered half-flux quanta at the 
corners of an unconnected zigzag structure (T=4.2K). The corner to corner distance is 
40 fim in both (a) and (b) . Reprinted figure with permission from Macmillan Publishers 
Ltd: H. Hilgcnkamp, Ariando, H.J.H. Smilde, D.H.A. Blank, G. Rijnders, H. Rogalla, 
J.R. Kirtley, and C.C. Tsuei, [Nature 422, 50 (2003)] Copyright 2003. 

shifts is shorter than the Josephson penetration depth, the total magnetic flux associated 
with each will be less than $o/2 [181] : a high density of facets can result in "splintered" 
Josephson vortices with fractions of $o of flux |182j . 

Facetted (zigzag) junctions can also be fabricated by design, using for example 
a ramp-edge YBCO/Nb junction technology |183j . An example is shown in Figure [T4"| 
which shows scanning SQUID microscope images of two facetted junctions |184j . In Fig. 
H4"a the corner junctions are connected by superconducting material. In this case there 
is a strong interaction between the half-flux quantum Josephson vortices upon cooling, 
so that they have a strong tendency to align with alternating signs of the circulating 
supercurrents and spontaneous magnetizations. In Fig. [T4"b the corner junctions have 
been isolated from each other by etching channels between them. In this case the half- 
flux quantum Josephson vortices all have the same alignment, even when cooled in a 
relatively small field. 

Lombardi et al. [185], and Smilde et al. [186] have mapped out the in-plane 
momentum dependence of the superconducting gap amplitude of YBCO by measuring 
the critical current as a function of junction normal angle for a series of biepitaxial 
YBCO grain boundary junctions, and YBCO/Nb ramp junctions respectively. The 
latter experiments, performed in part on untwinned epitaxial YBCO films, were able to 
determine the anisotropy of the gap between the a and b in-plane crystalline directions. 
However, these experiments were only sensitive to the amplitude of the superconducting 
order parameter, not its phase. Phase sensitive measurements of the momentum 
dependence of the in-plane gap in YBCO were performed by Kirtley et al. |187] by 
imaging the spontaneously generated magnetic flux in a series of YBCO-Nb two-junction 
rings, with one junction angle normal relative to the YBCO crystalline axes held fixed 
from ring to ring, while the other junction normal angle varied. Some results from this 
study are shown in Figure [Tj3 The rings either had n$o (integer flux quantization) or 
(n+ l/2)$o (half-integer flux quantization) of flux in them, n an integer. The transition 
between integer and half- integer flux quantization happened at junction angles slightly 
different from 45°+n x 90°, as would be expected for a pure d x 2_ y 2 superconductor, 
because of the difference in gap amplitudes between the crystalline a and b phases. In 
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Figure 15. (a) Schematic of YBCO-Nb two-junction rings fabricated to phase 
sensitively determine the momentum dependence of the in-plane gap in YBCO. (b) 
SQUID microscope image, taken with a square pickup loop 8 /jm on a side, of a square 
area 150 /zm on a side centered on one of the rings, cooled and imaged in zero field at 
4.2K. (c) The outer circle of images, taken after cooling in zero field, has a full-scale 
variation of 0.04 $o flux through the SQUID: the inner ring, taken after the sample 
was cooled in a field of 0.2 /iT, has a full scale variation of 0.09 $o- The rings cooled 
in zero field had either or $o/2 of flux in them; the rings cooled in 0.2 /iT had 
$0/2, $0 or 2$o of flux, as labelled. Reprinted figure with permission from Macmillan 
Publishing Ltd: J.R. Kirtley, C.C. Tsuei, Ariando, C.J.M. Verwijs, S. Harkema and 



H. Hilgenkamp, Nature Physics 2, 190 (2006) Copyright 2006. 



addition, careful integration of the total flux in the rings showed that any imaginary 
component to the order parameter, if present, must be small. 

3.1.4- Ring arrays Arrays of superconducting rings are of interest as a physical 
analogue of a spin: for conventional superconducting rings the states with shielding 
supercurrent flowing clockwise or counter-clockwise are degenerate at a magnetic flux 
bias of $o/2 (half of a flux quantum threading each ring). Two neighboring rings can 
interact magnetically, making them a physical analogue of the Ising spin. Frustration 
can be introduced into the system by choosing the arrangement of the rings: square and 
hexagonal lattices are geometrically unfrustrated, while triangular and Kagome lattices 
are frustrated. 

Davidovic et al. used scanning Hall bar microscopy to study two-dimensional arrays 
of closely spaced Nb rings (see Fig. [16]) |144[ 1145] . In these experiments Davidovic et al. 
found that it was possible to find small islands of antiferromagnetically ordered "spins" , 
but no perfect ordering, even in arrays with no geometrical frustration. One source of 
disorder in arrays of conventional superconducting rings is variation in area from ring to 
ring from lithographic variations. Conventional rings must be flux biased to near $o/2 
flux with high precision (see Fig. fT6l) . and area variations lead to flux variations. 

The same technology used for the pairing symmetry tests of Fig. [15] were used to 
form two-dimensional arrays of 7r-rings, superconducting rings with an intrinsic phase 
shift of 7r [184] . An example of SQUID microscope imaging of several 7r-ring arrays is 
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Figure 16. A sequence of field-cooled images of a honeycomb lattice of 1/xm diameter 
hexagonal Nb thin film rings, taken in increasing applied fluxes near an applied flux 
per ring of $o/2- The flux increases from left to right starting at the upper left corner 
with 0.4913$o and ending in the lower right corner at 0.5066<£>o- Reprinted figure 
with permission from D. Davidovic, S. Kumar, D.H. Reich, J. Siegel, S.B. Field, R.C. 
Tiberio, R. Hey, and K. Ploog, [Phys. Rev. B 55, 6518 (1997)] Copyright 1997 by the 
American Physical Society. 

shown in Figure [17] |188j . 7r-rings have two degenerate states at zero flux bias, therefore 
eliminating the need for precise flux biasing of the rings, and thereby eliminating 
slight differences in the sizes of the rings as a source of disorder [188J. Although tc- 
ring arrays showed more negative values of the bond-order a, corresponding to greater 
antiferromagnetic ordering, than arrays of conventional rings, perfect ordering was never 
observed in either conventional or 7r-ring arrays, possibly because of the large number 
of nearly degenerate states in this system |189] . 

3.1.5. Tests of the interlay er tunneling model One candidate mechanism for 
superconductivity at high temperatures in the cuprate perovskite superconductors is 
the interlayer tunneling model, in which the superconductivity results from an increased 
coupling between the layers in the superconducting state [1901 11911 1192] 1193] 1194] . An 
essential test of this theory is the strength of the interlayer Josephson tunneling in 
layered superconductors. For the interlayer tunneling model to succeed, the interlayer 
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Figure 17. SQUID microscopy images of four electrically disconnected arrays of 
7r-rings with 11.5 /im nearest neighbor distances. These images were taken at 4.2 
K with a 4 /im diameter pickup loop after cooling in nominally zero field at 1-10 
mK/s. The bond-order a is a measure of the anti-ferromagnetic order in the arrays. 
The geometrically unfrustrated square and honeycomb arrays could have perfect 
antiferromagnetic correlation, which would correspond to er — —1. The minimum 
possible bond-order for the frustrated Kagome and triangle lattices is a = —1/3. 
Reprinted figure with permissi on from J.R. Kirtley, C.C. T suei, Ariando, H.J.H. 
Smilde, and H. Hilgenkamp, |Phys. Rev. B 72, 214521 (2005)[ Copyright 2005 by the 
American Physical Society. 

coupling in the superconducting state must be sufficiently strong to account for the large 
condensation energy of the cuprate superconductors. The best materials for testing this 
requirement are Tl 2 Ba 2 Cu0 6+ 5 (Tl-2201) and HgBa 2 Cu0 4+ 5 (Hg-1201), which have high 
critical temperatures (T c ^90 K) and a single copper oxide plane per unit cell. Scanning 
SQUID microscope images were made of interlayer Josephson vortices emerging from 
the a - c face of Tl-2201 [TJS] (Fig. EES} and Hg-1201 USE]. The extent of these vortices 
in the a-axis direction is set by the interlayer penetration depth A c . Some spreading 
of the vortex fields occurs as the crystal surface is approached |197j . but not enough 
to affect the qualitative conclusions from these studies: the experimentally determined 
interlayer penetration depth A c is an order of magnitude larger than required by the 
interlayer tunneling model. These conclusions are in agreement with measurements of 
the Josephson plasma frequency in these materials |198j . 

3.1.6. Pearl vortices Vortices in thin superconductors satisfying d « Xl, where d is 
the thickness and Xl is the London penetration depth, first described by Pearl |199j . 
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Figure 18. A Scanning SQUID microscope image of an intcrlayer Josephson vortex 
emerging from the ac face of a single crystal of the single-layer cuprate superconductor 
Tl2Ba2Cu06+<5, imaged at 4.2K with a square pickup loop 8.2 ^m on a side. B A 
2-dimensional fit to the vortex in A with A c , the penetration depth perpendicular to 
the CuC>2 planes, equal to 18/im and zo, the spacing between the SQUID pickup loop 
and the sample surface, equal to 3.6 /im. C The difference between the data and 
the fit. D Cross sections perpendicular and parallel to the long axis of the vortex, 
offset for clarity. From K A. Moler, J.R. Kirtley, D.G. Hinks, T.W. Li, and Ming Xu, 
|Science 279, 1193 (1998)| Reprinted with permission from A A AS. 

have several interesting attributes. The field strengths H z perpendicular to the films 
diverge as 1/r at distances r << A, where A = 2\ 2 L /d is the Pearl length, in Pearl 
vortices, whereas in Abrikosov (bulk) vortices the fields diverge as ln(r/AL) [200J. Since 
in the Pearl vortex much of the vortex energy is associated with the fields outside of 
the superconductor, the interaction potential Vi nt (r) between Pearl vortices has a long 
range component Vi nt ~ A/r for r >> A |199j . unlike Abrikosov vortices, which have 
only short range interactions. Although it was originally believed that the interaction 
between Pearl vortices Vi n t ~ ln(A/r) for r << A leads to a Berezinskii-Kosterlitz- 
Thouless (BKT) transition which is cut off due to screening on a scale A Kogan 
|201j reports that the BKT transition could not happen in thin superconducting films 
of any size on insulating substrates because of boundary conditions at the film edges 
that turn the interaction into a near exponential decay. Grigorenko et ai, using SHM, 
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Figure 19. SQUID microscope image and cross-sectional data along the positions 
indicated by the dashed line in (a), (c) of Pearl vortices trapped in two artificially 
layered (Bao.gNdo.iCuO^+i)™/ (CaCuO^n superconducting films. The SQUID pickup 
loops were a square 7.5 /im on a side (a),(b), and an octagon 4/im on a side 
(c),(d) (schematics superimposed on image). The open symbols in (b),(d) are 
the cross-sectional data; the solid lines in (b), (d) are fits. Reprinted figure 
with permissi on from F. Tafu ri, J.R. Kirtley, P.G. Medaglia, P. Orgiani, and G. 
Balestrino, |Phys. Rev. Lett. 92, 157006 (2004)] Copyright 2004 by the American 
Physical Society. 



observed domains of commensurate vortex patterns near rational fractional matching 
fields of a periodic pinning array in a thin film superconductor, which could only form 
when the vortex- vortex interactions are long range [202] . 

Figure [19] shows scanning SQUID microscope images of Pearl vortices in ultra- 
thin [Bao.9Ndo.iCu0 2 +Jm/[CaCu0 2 ] n (CBCO) high-temperature superconductor films. 
Pearl vortices can be magnetically imaged in such thin films, although the Pearl length 
can be as large as 1 mm, because of the strong 1/r divergence of the Pearl vortex 
magnetic fields at the vortex core: The thin-film limit for the two-dimensional Fourier 
transform of the z-component of the field from an isolated vortex trapped in a thin film 
is given bv [T991I203] : 

*.<*.*) = ttS. (25) 



where z is the height above the film, k = + ky, A = 2\ ab /d, X a b is the in-plane 
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Figure 20. Scanning Hall probe microscope images of a 4 /im x 4/im square Pb 
film at applied fields of (a) 4.6 G, (c) 6.6 G, (e) 8.6 G, and (g) 10.6 G. [(b), (d), (f), 
and (h)] The results of 2D monopole fits to (a), (c), (e), and (g) respectively. The 
broken lines denote the shape of the square. The locations of the vortices are shown as 
white points. Reprinted figure with permission from T. Nishio, Q. Chen, W. Gillijns, 
K. De Keyser, K. Vervaeke, and V.V. Moshchalkov, |Phys. Rev. B 77, 012502~(2 008) 
Copyright 2008 by the American Physical Society. 



penetration depth, d is the film thickness and <3>o — h/2e. 

The open circles in Fig. HWb.d) are cross-sections through the data of Fig. [TWa.c). 
The solid lines are fits to the data using Eq. [251 using the Pearl length A as a fitting 
parameter. The in-plane London penetration depth X a b of these thin films can be inferred 
from A if d is known. 

3.1.7. Mesoscopic superconductors Vortices in superconductors with dimensions 
comparable to the superconducting coherence length £ have been predicted to nucleate 
spontaneously in a number of interesting configurations, such as shells [204] . giant 
vortices |205j . and anti- vortices, so as to preserve the symmetry of the sample 
[206[ 1207] . Vortex trapping patterns in mesoscopic superconducting samples have been 
studied using Bitter decoration |208] . scanning SQUID microscopy [209J, ballistic Hall 
magnetometry If) .1 . transport |210] . and Hall bar microscopy. An example of Hall 
bar microscopy results is shown in Figure [20j Nishio et. al. [211] report results in 
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Figure 21. Scanning SQUID microscope images of flux domains in Sr2Ru04 at 
T=0.36 K after field cooling at various fields. In all cases, the magnetic field amplitude 
applied along the c axis (Hj_) was kept constant at 2 G while the in-plane field (H a b) 
was (a) G, (b) 5 G, (c) 10 G, (d) 50 G. The imaging area is 31 ^mx 17 urn. The 
field scale in G is shown on the right; dark regions are superconducting vortex-free 
regions. Reprinted figure with permission from V.O. D olocan, C. Veauvy, F. Servant, 
P. Lejay, K. Hasselbach. Y. Liu, and D. Mailly, |Phys. Rev. Lett. 95, 097004 (2005)| 
Copyright 2005 by the American Physical Society. 

agreement with Ginzburg-Landau calculations for temperatures considerably lower than 
the superconducting transition temperature. 

3.1.8. Sr 2 RuOi Sr 2 Ru04, with the same crystal structure as L^CuO/i, the parent 
compound for the first class of cuprate perovskite compounds shown to exhibit high- 
temperature superconductivity [2J, was discovered in 1994 to be superconducting at 
about IK [5J. It is thought to be a p x ± ip y -w&ve, triplet pairing superconductor with 
broken time-reversal symmetry. One of the consequences of this pairing symmetry 
state is spontaneously generated supercurrents, and consequent large magnetic fields 
at surfaces and boundaries between p x + ip y and p x — ip y domains |212] . Although 
there is evidence for characteristic magnetic fields of 0.5G in Sr2Ru04 from muon spin 
resonance experiments |213j . and evidence for broken time-reversal symmetry at the 
surface of Sr 2 Ru0 4 has been obtained using Sagnac interferometry [TS], no evidence 
for spontaneously generated magnetic fields at the surfaces of Sr2Ru04 samples has 
been found in scanning Hall bar and scanning SQUID experiments, with much higher 
sensitivity than the magnetic fields originally predicted |214[ I215[ 1216] . A possible 
explanation for this failure to observe spontaneous magnetization directly is that the 
p x ± ip y domains are small, so that the fields from closely spaced domain boundaries 
nearly cancel each other. The experimental information on p x ± ip y domain sizes is 
not consistent, with the first phase sensitive pairing symmetry experiments on Sr2Ru04 
being consistent with domain sizes of order 1 mm |217j . while later phase-sensitive 
results are consistent with domain sizes 1/im or less |218j . and the Sagnac interferometry 
being consistent with domain sizes intermediate between these two sizes [18]. Recently 
Raghu et al. have suggested that the superconductivity in Sr 2 Ru04 arises primarily in 
quasi-one- dimensional bands, and is therefore not expected to generate experimentally 
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Figure 22. Scanning Hall bar image of a single crystal sample of ErN^I^C. 
After cooling below T c in a weak field [2.4 Oe in (a), (b), (d), (e), 18 Oe in (c)], 
a vortex distribution consistent with random pinning is observed on the ab face 
[(a) 7.3 K]. Upon reducing the temperature below Tjv, the vortices spontaneously 
organize along twin domain walls along the [110] [(b) 5.3 K] or [110] [(c) 4.2 K] 
direction. The pattern gradually disappears as the temperature is raised again 
[(d) 5.7 K, (e) 6.0 K, different cycle than (a) and (b)]. Reprinted figure with 
permissio n from H . Bluhm , S.E. S ebastian, J.W. Guikema, I.R. Fisher, and K.A. 
Moler, |Phys. Rev. B 73, 014514 (2006)| Copyright 2006 by the American Physical 
Society. 

detectable edge currents |219j . 

Dolocan et al. [2201 1221] have magnetically imaged single crystals of Sr2Ru04 
using a scanning /xSQUID mounted on a tuning fork in a dilution refrigerator. This 
allowed combined magnetic and topographic information on a length scale of 1 micron, 
at temperatures below 400 mK. They found coalescence of vortices and the formation 
of flux domains (see Fig. |2TI) . The formation of lines of vortices in a field tilted relative 
to the crystalline c-axis [222] can be attributed in this case to the large anisotropy in 
Sr2Ru04, resulting in an attractive interaction between vortices in a plane defined by 
the c-axis and the magnetic field direction |221j . The vortex coalescence observed when 
the applied magnetic field is parallel to the c axis may be due to weak intrinsic pinning 
at domain walls, and the existence of a mechanism for bringing vortices together that 
overcomes the conventional repulsive vortex-vortex interaction. 

3.1.9. Magnetic superconductors Magnetic order and superconductivity are competing 
orders since the Meissner state excludes a magnetic field from the bulk and 
superconductivity is destroyed at sufficiently high fields. However, superconductivity 
and magnetism can coexist [2231 12241 12251 1226] if the orientation of the local magnetic 
moments varies on a length scale shorter than the superconducting penetration depth A, 
or if the field generated by the magnetization is carried by a so-called spontaneous vortex 
lattice [227J. The former order has been observed in ErRh 4 B 4 and HoMo 6 S 8 |228[ 1229]. 
and indirect evidence for the latter has been reported in UCoGe [230]. Several Hall 
bar studies have been performed on artificial ferromagnetic-superconducting hybrid 
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Figure 23. Cross-sections and images of single (bottom right image, blue cross- 
section) and double (middle right image, red cross section) flux quanta trapped in a 
50/im x 50^m square of single-crystal YBa2Cu3 06.35 (T C =6.0K) with a 10 /zm hole 
drilled using a focussed ion beam (top right image). When the sample was heated to 
5.6K for 1 sec and re-cooled to 2K, these flux quanta escaped with no sign of the vortex 
memory associated with visons. Reprinted figure with permission from Macmillan 
Publishing Ltd: D.A. Bonn, J.C. Wynn, B.W. Gardner, Y- J. Lin, R. Liang, W.N. 
Hardy, J.R. Kirtley and K.A. Moler, [Nature 414, 887 (2001)| Copyright 2001. 



structures [2311 12321 12331 1234] . Fig. [22] [235] shows scanning Hall bar microscopy 
measurements of ErN^E^C, which has a superconducting T c of ~ 11 K, becomes 
antiferromagnetic below TV ~ 6 K, and exhibits weak ferromagnetism below T^wfm ~ 
2.3 K [236J. Bluhm et al. [235J found that in this material the superconducting 
vortices spontaneously rearranged to pin on twin boundaries upon cooling through 
the antiferromagnetic transition temperature, and that a weak random magnetic signal 
appears in the ferromagnetic phase. 

3.2. Fluxoid dynamics 

3.2.1. Limits on spin-charge separation One explanation for the peculiar normal state 
properties and high superconducting transition temperatures of the high-T c cuprate 
perovskites is that there exists a new state of matter in which the elementary excitations 
are not electron-like, as in conventional metals, but rather the carriers "fractionalize" 
into "spinons", chargeless spin 1/2 fermions, and "chargons", bosons with charge +e 
[237J. Senthil and Fisher proposed a model for such a separation in two dimensions 
with sharp experimental tests [238] : in conventional superconductors a Cooper pair 
with charge 2e circling around a vortex with h/2e of total magnetic flux experiences 
a phase shift of 2ir. However, a chargon would pick up only a phase shift of ir. In 
order for the chargon wave function to remain single-valued, it must be accompanied 
by a "vison" that provides an additional 7r phase shift. In the Senthil-Fisher test, a 
superconducting ring with flux trapped in it is warmed through the superconducting 
transition. If the number of vortices in the ring is odd, then it will also contain a vison 
that could persist above T c . If the ring is then recooled in zero field before the vison 
escapes, it will spontaneously generate a single flux quantum with random sign. On 
the other hand, if the number of vortices in the ring is even, no such "vortex memory" 
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Figure 24. SQUID microscope images of a 12 x 12 array of 20 fim inside 
diameter, 30 /mi outside diameter thin film rings of MoaSi, cooled in zero field 
through the superconducting T c at various rates. The resulting ring vortex 
numbers were consistent with thermally activated spontaneous fluxoid formation. 
Reprinted figure with permission from J.R. Kirtley, C.C. Tsuei, and F. Tafuri, 
|Phys. Rev. Lett. 90, 257001 (2003)"] Copyright 2003 by the American Physical 
Society. 

effect would exist. Bonn et al. |239] (Fig. [23]) tested for this effect in highly underdoped 
single crystals of YL^CuaOe+z- These crystals have a very sharp but low transition 
temperature favorable for such a test. They found no evidence for a vortex memory 
effect, and were able to put stringent upper limits on the energy associated with a 
vison. Some models in which superconductivity results from spin-charge separation 
predicts h/e fluxoids in materials with low superfluid density [24QJ I241[ 1242] . Wynn 
et al. |243j, using scanning SQUID and Hall bar microscopies, found no evidence for 
h/e vortices in strongly underdoped YBa 2 Cu30 6+a; crystals, also placing limits on spin- 
charge separation in these materials. 

3.2.2. Superconducting rings As we have seen in previous sections, the dynamics of 
fluxoids in a ring configuration are of interest as a model for the Ising spin system 
|144[ I145[ 1188] I3.1.4[ and for placing limits on spin-charge separation in underdoped 
cuprate superconductors [239J 13.2.11 Kirtley et al. studied fluxoid dynamics in 
photolithographically patterned thin film rings of the underdoped high-temperature 
superconductor Bi2Sr2CaCu20 8 +,5 using a scanning SQUID microscope, and concluded 
that their results could be understood in terms of thermally activated nucleation of a 
Pearl vortex in, and transport of the Pearl vortex across, the ring wall [244J. 
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Fluxoid dynamics in superconducting rings may also provide clues to physics on an 
entirely different scale: the early development of the universe. Immediately after the Big 
Bang a fundamental symmetry related the particles that mediated the electromagnetic, 
strong, and weak forces. However, as the universe cooled down this symmetry was 
broken in such a way that, for example, the W and Z bosons, which mediate the 
weak interaction, have mass, while photons, which mediate the electromagnetic force, 
do not. Kibble proposed an intuitive picture of this symmetry breaking of the early 
universe [245[ [246J, with the underlying idea that causality governed the number of 
defects because a finite time delay is needed for information to be transferred between 
different regions of a system. Zurek pointed out that Kibble's ideas could be tested 
by studying vortices in superfluids and superconductors |247] . As a superconductor or 
superfluid cools through its transition temperature different regions can nucleate into 
a state with the same order parameter amplitude ip but different phases 0. This will 
produce "kinks" or "defects" where would eventually (at low temperature) change by 
2tt. These kinks can be removed by annihilation with anti-kinks if the coherence length 
£ is long enough and the cooling rate is sufficiently slow. However, as the system cools 
this annihilation process proceeds less rapidly, until at some point the kinks "freeze-in" 
because part of the system cannot communicate with other parts quickly enough. In 
a superconducting ring a phase shift of 2n around the ring corresponds to a fluxoid 
trapped in the ring, so that the process of forming and freezing in kinks is also called 
spontaneous fluxoid generation. Kibble and Zurek showed that causality implies that the 
probability P of finding a spontaneous fluxoid in the ring (when cooled in zero external 
field) depends on the cooling rate Tq 1 as P oc {tq/tq)~ (J , where r is a characteristic 
time and a depends on how the correlation length £ and the relaxation time r vary 
with temperature close to the transition. For a superconducting ring and assuming a 
mean-field approximation a = 1/4. 

Tests of the Kibble-Zurek prediction have been made by looking for vortices in 
superfluid *He [2HJ and 3 He (233 I25T] . in superconducting films E5], and 
superconducting rings interrupted by Josephson junctions [253[ I254[ [255J. However, 
verification of the ideas of Kibble and Zurek using arguably the simplest system, 
superconducting rings, were not made until recently. Kirtley, Tsuei and Tafuri |256] 
(Fig. 124")) used a SQUID microscope to image an array of rings after repeated cooling 
in various magnetic fields and cooling rates to determine the probability of spontaneous 
fluxoid generation. They showed that a second mechanism prevailed, in which the final 
density of fluxoids depended on a balance between thermal generation and the relaxation 
rates of fluxoids |257[ I258[ 1259] . They argued that the thermal fluctuation mechanism 
is complementary to the "causal" mechanism and should be considered in attempts to 
understand phase transitions, both in the laboratory and in the early universe. The 
Kirtley et al. experiments |256] were done on rings with parameters which favored 
the thermal activation mechanism over the Kibble-Zurek (causal) mechanism. Recently 
Monaco et al. [260] have studied superconducting rings with parameters more favorable 
to the causal mechanism, using a geometry in which the flux state of a single ring can 
be manipulated and sensed rapidly, and find results consistent with the Kibble-Zurek 
prediction, taking into account the fact that the ring's circumference was much smaller 
than the coherence length at the temperature at which the fluctuations were frozen in. 

3.2.3. High-T c grain boundaries Grain boundaries in the cuprate high-T c supercon- 
ductors are widely used for devices and fundamental studies, and govern the transport 
properties of superconductors with technological applications |178j . It is therefore im- 
portant to understand the mechanism of dissipation in transport across grain bound- 
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Figure 25. Scanning Hall bar measurements of vortex noise in a 24° symmetric 
YBCO grain boundary (GB). (a) Sketch of the sample. The bridge width is 50 /im. 
(b,c) Magnetic image of the GB area (b) at 10 K and applied current I = 29 mA 
and (c) at 10 K after thermal cycling to 70 K with 1 = 0, showing individual vortices 
(red) and antivortices (blue). The regions imaged are outlined in orange in (a). (d,e) 
The Hall probe signals vs. time at 10 K for different currents applied to the GB 
measured at two positions marked in Fig. lb. The switching events are vortices 
passing under the probe. Red circles indicate three-level switching events. Reprinted 
figure with permission from B. Kalisky, J.R. Kirtley, E.A. Nowadnick, R.B. Dinner, 
E. Zeldov, Ariando, S. Wenderich, H. Hilgenkamp, D.M. Feldmann, and K.A. Moler, 
|Appl. Phys. Lett 94, 202504 (2009)] Copyright 2009 by the American Institute of 
Physics. 



aries. Transport across grain boundaries is usually studied with a voltage threshold 
on the order of nanovolts. At these voltages millions of vortices traverse the bound- 
ary per second. Much smaller voltages can be measured by magnetically imaging the 
vortices in the grain boundaries. Kalisky et al. [261] performed such measurements on 
grain boundaries produced by epitaxial growth of the cuprate high-T c superconductor 
YBa 2 Cu307_5 on bicrystals of SrTiC>3 using a large scanning area Hall bar microscope 
|84j . They observed (Figure l25l) telegraph noise in the Hall bar signal when the sensor 
was directly above the grain boundary, which they attributed to the motion of vortices. 
They inferred the voltage across the grain boundary from the frequency of the telegraph 
noise, with voltage sensitivity as small as 2xlO~ 16 V. The dependence of the inferred 
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Figure 26. SQUID signal (left axis) and ring current (right axis) as a function of 
applied flux $ a for two aluminum rings, both with thickness d = 60 nm and annulus 
width w — 110 nm. The fluctuation theory (dashed red) was fit to the data (blue) 
through the dependence of / r i ng on $ a and temperature. (A to C) ring radius R = 
0.35 /im, fitted T c ($ a =0) = 1.247 K, and 7 = 0.075. The parameter 7 characterizes 
the size of the ring. The green line is the theoretical mean field response for T = 
1.22 K and shows the characteristic Little-Parks line shape, in which the ring is not 
superconducting near $ a = $o/2. The excess persistent current in this region indicates 
the large fluctuations in the Little-Parks regime. (D) R = 2 /im, fitted T c ($ a =0) — 
1.252 K, and 7 = 13. F rom N.C. Koshnick, H. Bluhm, M.E. Hubcr, and K.A. Moler, 
|Science 318, 1440 (2007)| Reprinted with permission from A A AS. 



grain boundary voltage on current followed either an exponential V oc e ml l Ic or a power 
law V oc I n dependence with high values of m or n. These results were qualitatively dif- 
ferent from grain boundary transport measurements at higher voltages |178[ I262[ |263j , 
and could not be fully explained using existing models. 

3.3. Local susceptibility measurements 

3.3.1. Superconducting fluctuations Experimental knowledge of superconducting 
fluctuations in one dimension has largely been derived from transport measurements 
[264J, which require electrical contacts and an externally applied current. The SQUID 
susceptometer, along with a ring geometry for the superconductor, allows the study 
of fluctuation effects, without contacts, in isolated, quasi one-dimensional rings in the 
temperature range where the circumference is comparable to the temperature-dependent 
Ginzburg-Landau coherence length £(T). In such measurements, an example of which 
is shown in Fig. [26] [265] , a magnetic field is applied to the ring using a single-turn field 
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coil, co-planar and concentric with the pickup loop, both of which are integrated into 
the SQUID sensor [72]. The response of the ring to this applied flux is measured by the 
pickup loop. The SQUID sensor is in a carefully balanced gradiometer configuration, 
such that it is insensitive to the fields applied by the field coil [73] • Nevertheless, further 
background subtraction, by comparing the mutual inductance between the field coil 
and SQUID with the pickup loop close to the ring with that with the pickup loop at a 
distance, is required to measure the response field due to the ring, which can be 7 orders 
of magnitude smaller than the applied field. A scanned sensor allows multiple rings to 
be measured in a single cooldown. The experimental results agree with a full numeric 
solution of thermal fluctuations in a Ginzburg-Landau framework that includes non- 
Gaussian effects |266[ 1267] for all of the rings for which calculations were numerically 
tractable. This is in contrast to previous work on a single Al ring [268], which disagreed 
strongly with theory. 

Similar techniques have allowed the measurement of spontaneous persistent currents 
in normal metal rings [269]. 

3.3.2. Stripes The pnictide class of superconductors [8] have the highest critical 
temperatures (57K) observed for a non-cuprate superconductor, multiband Fermi 
surfaces, and at certain doping levels a paramagnetic to antiferromagnetic as well 
as a tetragonal to orthorhombic transition above the superconducting transition 
temperature. In Ba(Fei_ a; Co; C )2As2 (Co doped Ba-122) and other members of the 
122 family (AFe2As2 with A=Ca, Sr, and Ba), doping causes the spin-density- wave 
transition temperature and the structural transition temperature to decrease |270[ 1271] 
falling to zero at or near the doping where the highest T c occurs, suggesting the 
importance of lattice changes in determining transport properties. Both experiment 
[272[ 12731 1274j and theory [275] have suggested a close relationship between structural 
and magnetic properties, leading some authors to describe the lattice and spin density- 
wave transition by a single order parameter [276[ 1277] . In addition, structural strain 
appears to play a significant role in the superconductivity. For example, in the 122 
compounds, small amounts of non- hydrostatic pressure can induce superconductivity 
[2781 12791 [2501 [28H [282]. In addition, there is evidence that the structural perfection of 
the Fe-As tetrahedron is important for the high critical temperatures observed in the 
Fe pnictides [28311254] . 

Kalisky et al. [285] have shown that Co doped Ba-122 has stripes of enhanced 
susceptibility below the superconducting transition temperature (see Fig. [27]). These 
stripes are resolution limited using a SQUID pickup loop with an effective diameter 
of 4/zm, and are believed to be associated with boundaries between twins- crystallites 
with their in-plane a, b axes rotated by 90°. The amplitude of the susceptibility stripes 
becomes larger as the superconducting transition temperature is approached from below. 
Since the susceptibility signal for a homogeneous superconductor becomes larger as 
the penetration depth becomes shorter and the superfluid density becomes larger, it 
seems reasonable to associate the susceptibility stripes with shorter penetration depths 
and enhanced superfluid densities on the twin boundaries. Twin boundaries have been 
associated with enhanced superfluid density in conventional superconductors [286j . but 
the influence of twin planes on superconductivity in, for example, the cuprates is less 
clear [2871 [2881 [2891 121M I29T] . 

It is difficult to calculate analytically the effect of a plane of reduced penetration 
depth (enhanced superfluid density) on the susceptibility observed in scanning SQUID 
measurements. Kirtley et al. [292J have done a finite-element calculation in the 
appropriate geometry. Because of the uncertainty in the width of the region with 
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Figure 27. Local susceptibility image in underdoped Ba(Fei_ a; Co; C )2As2, indicating 
increased diamagnetic shielding on twin boundaries. (a) Local diamagnetic 
susceptibility, at T=17 K, of the ab face of a sample with x=0.051 and T c =18.25 
K, showing stripes of enhanced diamagnetic response (white). In addition there is a 
mottled background associated with local T c variations that becomes more pronounced 
as T — > T c . Overlay: sketch of the scanning SQUID sensor. The size of the pickup loop 
sets the spatial resolution of the susceptibility images, (b) and (c) Images of the same 
region at (b) T=17.5 K and (c) at T=18.5 K show that the stripes disappear above 
T c . A topographic feature (scratch) appears in (b) and (c). (d) Crystal structure of 
a unit cell of the parent compound BaFe2As2. (e) Top view of the FeAs layer with 
tetragonal symmetry, and (f) an exaggerated view of the orthorhombic distortion that 
occurs at low temperatures, (g) Schematic of a possible arrangement of spins across a 
twin boundary in the anti-ferromagnetic state. Reprinted figure with permission from 
B. Kalisk y, J.R. Kirtley, J.G. Analytis, J .-H. Chu, A. Vailionis, I.R. Fisher, and K.A. 
Moler, |Phys. Rev. B 81, 184513 (2010)| Copyright 2010 by the American Physical 
Society. 

enhanced superfluid density, it is also difficult to estimate the size of the enhancement in 
the Cooper pair density from the scanning susceptibility measurements. Nevertheless, 
Kirtley et al. estimate an enhancement in the two-dimensional Cooper pair density in 
the range between 10 19 and 10 20 m~ 2 . For comparison, the two-dimensional electron 
liquid |293] at the LaA103-SiTi03 interface, which exhibits superconductivity at about 
0.2K |294] . has a carrier concentration of about 10 17 m -2 . The temperature dependence 
of the stripe amplitude can be best fit by assuming that the twin boundary has a different 
critical temperature than the bulk, although no stripes were observed above the bulk 
T 

Kalisky et al. also observed that vortices tended to avoid pinning on the twin 
boundaries, and that they were difficult to drag over the twin boundaries using either a 
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SQUID susceptometer sensor or a magnetic force microscope tip [292J. The enhancement 
of supercurrent density, as well as the barrier to vortex motion presented by the twin 
boundaries, provide two mechanisms for enhanced critical currents in twinned samples. 
Prozorov et al. |295j have noted an enhancement of the critical current of slightly 
underdoped single crystals of Ba(Fei_ x Co x ) 2 As 2 , which they associate with twinning. 



3.3.3. Spinlike susceptibility There appears to be a component of 1/ f a (a ~ 1) noise 
in superconducting devices, which is apparently due to magnetic flux noise, and is 
remarkably universal [2961 1297] . Recently Koch, DiVincenzo, and Clarke (KDC) have 
proposed that this noise is due to unpaired, metastable spin states [44J. The proposed 
spin states are bistable, with a broad distribution of activation energies for changing 
their spin orientation. When the spin state changes it changes the amount of magnetic 
flux coupling into the superconducting device. Each contributes a Lorentzian with a 
particular cutoff frequency to the total noise; adding up a large number of Lorentzians 
with a broad distribution of cutoff frequencies results in noise with 1/f distribution [43J. 
KDC estimate that a density of unpaired spin states of about 5xl0 17 m~ 2 is required to 
explain universal 1/f noise. 

Bluhm et al. [298J have found a signal in scanning SQUID susceptometry 
measurements at low temperatures in a number of materials, both metals and insulators, 
even including Au, that has a paramagnetic response with a temperature dependence 
consistent with unpaired spins. This susceptibility has a component that is out of phase 
with the applied field, implying that it could contribute to 1 / /-like magnetic noise. The 
density of these spin states is estimated to be in the range of 10 17 m -2 , consistent with 
the KDC estimate. This implies that scanning SQUID susceptometry could be used as 
a diagnostic for determining which materials and processes contribute most strongly to 
1/f noise in superconducting devices. 



3.4- Penetration depths 

Penetration depths have been inferred from fitting the magnetic images of 
superconducting vortices for a number of years [79j [801 I299[ 1321 13001 I301[ 1302] . An 
alternate means is to measure the mutual inductance between the SQUID pickup loop 
and a field coil integrated into the sensor in SSM [72] (Fig. 1281) . The 2-dimensional 
Fourier transform of the z-component of the response field in the pickup loop, with a 
current J in a circular ring of radius R oriented parallel to, and a height z above a 
homogeneous bulk superconductor with penetration depth A is given by [303J 

B z (k) = -tt^IR J x {kR)e- 2kz , (26) 



where q = yk 2 + l/X 2 , and J\ is the order 1 Bessel function of the first kind. Similarly 
the response field for a thin film superconducting sample with Pearl length A is given 
by [303] 

— 2kz 

B z {k) = -TTfioIRJ^kR)^——. (27) 

1 + kA 

To a good approximation, if z » A the magnetic fields resulting from the screening of 
the field coil fields by the sample act as if they are due to an image coil spaced by a 
distance 2h e ff = 2(z + A) from the real coil, where z is the physical spacing between 
the sample surface and the field coil. A SQUID sensor with a single circular field coil 
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Figure 28. (a) Geometry of a SQUID susceptometer. (b) Mutual inductance between 
the field coil and the pickup loop, as a function of the spacing z between the SQUID 
substrate and a thin CBCO film. The symbols are data, taken with various alternating 
currents through the field coil. The solid line is modeling using Eq. [331 with A=11.8/im. 
(b) Comparison of the Pearl length A for a number of CBCO samples using fitting of 
SQUID magnetometry images of vortices (vertical axis) vs susceptibility measurements 
(horizontal axis). Reprinted figure with permission from F. Ta furi, J.R. Kirtley, 
P.G. Medaglia, P. Orgiani, and G. Balestrino, |Phys. Rev. Lett. 92, 157006 (2004)] 
Copyright 2004 by the American Physical Society. 



of radius R and a co-planar, coaxial circular pickup loop of radius a oriented parallel to 
a homogeneous bulk sample surface will have a mutual inductance [304J 

M = (i - wrk^) ' (28) 

The solid line in Fig. 128T b) is obtained by numerically integrating the 2-D Fourier 
transform of Eq. ([27)1 over the area of the pickup loop, for various values of z, and 
fit to the data by varying A. Figure 1281(c) compares the values obtained for the Pearl 
lengths for a number of the CBCO samples using magnetometry and susceptometry 
methods. The two methods agree within experimental error over the range of Pearl 
lengths present. Since the fitting to the Pearl vortex images was done assuming each 
vortex has $o of total flux threading through it, rather than a fractional value |3U5] . this 
agreement means that the superconducting layers are sufficiently strongly Josephson- 
coupled that it is energetically favorable for the vortex flux to thread vertically through 
the superconducting layers, as opposed to escaping between the layers. 

Figure [29] shows results from a study using scanning SQUID susceptometry of 
the penetration depth of the pnictide superconductor LaFePO [304J. This study 
used the technique of measuring the mutual inductance between the SQUID field coil 
and pickup loop as a function of spacing z, and then repeating the measurement at 
fixed z while varying the temperature. Since the mutual inductance is a function of 
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Figure 29. (a,b) Optical micrographs of two single crystals of the pnictide 
superconductor LaFePO. (c) Susceptibility scan of #1 at T=0.4K. (d) Change in 
hefi = h + A between 0.4 and 3K over sample #2. (e,f ) Maps of local T c over the same 
areas as in (c) and (d). The crosses indicate points where AA(T) data were collected. 

(g) AA of samples #1, and #2. The dashed lines are fits between 0.7 K < T < 1.6 K. 

(h) Black lines are possible superfluid densities for LaFePO sample #1, point 2, with 
different A(0). Shaded area: superfluid density of YBa2Cu306.gg from [306] . [307] . The 
width of the shaded area reflects uncertainty in A(0). Reprinted figure with permission 
from C.W. Hicks, T.M. Lipp man, M.E . Huber , J.G. Analytis, J.-H. Chu, A.S. Erickson, 



I.R. Fisher, and K.A. Moler, |Phys. Rev. Lett. 10 3, 127003 (2009) Copyright 2009 by 
the American Physical Society. 



z + A, the temperature dependence of A can be inferred from such measurements, even 
without a detailed model for the mutual inductance. Spatially resolved susceptibility 
measurements (Fig. I29b-f) showed that the effective height depended on the topography 
of the sample, with scratches, bumps, and pits effecting the measurements strongly. 
However, the measured penetration depth and the temperature dependence of the 
penetration depth was reproducible for positions more than 10/xm from surface 
irregularities. In this case the penetration depth had a linear temperature dependence at 
low temperatures, similar to the cuprate high temperature superconductors, indicative 
of well formed nodes in the energy gap. It should be noted that in this geometry it is 
X a b, the in-plane penetration depth, that controls the mutual inductance. 

Perhaps the most promising application of MFM to the study of superconductivity 
is for high spatial resolution, absolute measurements of penetration depths. There are 
two ways of inferring superconducting penetration depths from MFM measurements. 
The first is, as mentioned above, to fit MFM vortex images to a model with the 
penetration depth as a fitting parameter [30 11 1302] . While Nazaretski et al. modeled 
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Figure 30. Normalized superfluid density p s (T)/p s (0) = X ab (0) 2 / \ ab {T) 2 vs T for 
the pnictide superconductor Ba(Feo.g5Coo.o5)2As2. AA a b(T) is determined from MFM 
(squares) and SSS (diamonds) experiments by measuring the change in the diamagnetic 
response at fixed height. These values are offset to match the absolute value of \ a b(T) 
by fitting the MFM data to a truncated cone model (circles). The green solid line 
shows a fit to a two-band s-wave model. The width of the dashed band reflects the 
uncertainty in A a b(0). Inset: AA Q b(T) vs T at low T. Black dashed line: one gap s-wave 
model with <z=1.5 and A = 1.95T C . Magenta dashed line: AA(T) = cT 22 (c=0.14 
nm/K 2 - 2 ). Reprinted figure with permission from L. Luan, O.M. Auslaender, T.M. 
Lippman, C.W. Hicks, B. Kalisky, J.-H. Chu, J.G. Analytis, I.R. Fisher, J.R. Kirtlcy, 
and K.A. Moler, |Phys. Rev. B 81, 100501 (2010)| Copyright 2010 by the American 
Physical Society. 



their tip numerically [[301], Luan et al. used an analytical approach which is easy to 
fit to experimental data and not strongly affected by differences between tips [302J. 
A second method is to consider the interaction between the magnetic material in the 
tip and its image in the superconducting material due to Meissner screening of the tip 
fields. Following Xu et al. |3U8j . Lu et al. |3U9] wrote the force due to the interaction 
between the magnetic tip, approximated as a magnetic dipole, and its image in the 
superconductor as 



where ttlq is the tip's magnetic moment and z is the physical distance between tip 
and sample. They fit their measured MFM force derivative curves to the z » A 
approximation 




(29) 





(30) 



dz 



8n[z + A(T)] 5 ' 
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to obtain the temperature dependence of the penetration depth A for a single crystal of 
YBCO. Luan et al. [302] used a truncated cone tip model to write 

dF g (z,T) dF z (z,T) 



dz dz 

\z + X ab (T)^ [z + \ ab {T)f^2[z + \ ab {T)f\ 1 ] 

for the difference between the force derivative at a given height z and that at large 
distances, where A, which depends on the tip shape and coating, is determined from 
fitting experimental data at low temperatures, \ a b{T) is the in-plane component of the 
penetration depth, and ho is the truncation height of the tip. Model independent values 
for AA(T) were obtained by mapping a change in A to a change in z since they are 
equivalent. Absolute values for A(T) were obtained by fitting to Eq. EU Figure I3"U1 
shows experimental results for the superfluid density p s of the pnictide superconductor 
Ba(Feo.95Coo.o5)2As2 obtained by MFM and SSM measurements. The AA(T) results 
from MFM measurements were later confirmed by a bulk technique on similar samples 
[310J. The advantage of MFM is that both the absolute value and high sensitivity 
relative values of A can be obtained in the same measurement, so that superfluid density 
is obtained over the full temperature range and from it the gap structure of the order 
parameter. Luan et al. found a two full-gap model model fit the data well, consistent 
with the s +_ order parameter most frequently discussed in the pnictides. Luan et al. 
also found that p s was uniform on the submicron scale despite highly disordered vortex 
pinning. 



3. 5. Manipulation of individual vortices 

While there is a vast literature on the properties of superconducting vortices [111] . 
there has been relatively little work on manipulation of individual vortices by scanning 
magnetic probes. Such manipulation, performed with MFM |112U119l II 16113111 1312] I3l3] 
or SSM [314] . can directly measure the interaction of a moving vortex with the local 
disorder potential [312] . Keay et al. used MFM to observe sequential vortex hopping 
between sites in an array of artificial pinning centers [313] . An example of another such 
MFM study is shown in Figure [31] [312] . Here vortices in a high quality, detwinned 
sample of YL^CuaOg.ggi are fixed at low temperatures and high scan heights, but at 
higher temperatures and lower scan heights the tops of the vortices can be dragged 
by an attractive magnetic interaction with the tip. The vortices are dragged along 
the fast scan direction until a combination of pinning and elastic forces overcomes the 
attractive tip- vortex force. The displacement of the vortex in the slow scan direction is 
much longer than in the fast scan direction, demonstrating a "vortex wiggling" effect: 
alternating transverse forces markedly enhance vortex dragging. Vortex motion in both 
the slow and fast scan directions is stochastic. From such measurements Auslaender et 
al. [312J were able to infer the maximum lateral dragging force before the top of the 
vortices were depinned. They also found a marked in-plane anisotropy in the distance 
w the top of the vortex was dragged by the tip in the fast scan direction, larger than 
could be accounted for by the in-plane anisotropy of the vortex core and point pinning. 
They speculated that a likely source for extra pinning anisotropy is nanoscale clustering 
of oxygen vacancies along the Cu-0 chains in YBCO. 



44 




Figure 31. MFM imaging and manipulation of individual vortices in YBCO at T = 
22. 3K. a,b Schematic drawings of an MFM tip (triangle) that attracts a vortex (thick 
lines) in a sample with randomly distributed pinning sites (dots): at large heights a 
the force between tip and vortex is too weak to move the vortex; at small heights b 
the vortex moves right and then left as the tip rasters over it. c MFM image taken 
at z—420 nm (maximum applied lateral force F^ x m 6 nN). d z=170 nm (F^* x w 
12 nm). Inset: Scan taken at a comparable height at T=5.2K. e Line cut through the 
data in c along the dashed line, showing the signal from a stationary vortex (blue). 
Overlapping it is a line cut from the reverse scan (green), f Line cut through the data 
in d along the dashed line, showing a typical signal from a dragged vortex. Reprinted 
by permission from Macmillan Publishers Ltd: OM. Auslaender, L. Luan, E.W.J. 
Straver, J.E. Hoff man, N.C . Koshnick, E . Zeldov , D.A. Bonn, R. Liang, W.N. Hardy, 
and K.A. Moler, |Nature Physics 5, 35-39 (2008) | Copyright 2008. 



4. Conclusions 

In summary, there are now many types of scanning magnetic microscopies, each 
with their own strengths and weaknesses. Of the three covered here, SSM is the 
most sensitive for low resolution applications, but requires a cooled sensor. MFM 
has the highest spatial resolution, but also applies the largest fields to the sample. 
SHM can be used over a broader temperature range than SSM. A major thrust in 
development has been towards higher spatial resolution, which requires smaller sensors 
scanned closer to the sample. The scaling of resolution and sensitivity with sensor 
size and spacing depends both on the technique and the source of field. Many of 
the applications discussed here have involved the imaging of magnetic flux in special 
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geometries, unconventional superconductors, or both. Others have involved time 
dependence, through e.g. multiple cooldowns, cooldowns at varying cooling rates, or 
noise measurements. Some exciting new developments have been scanning susceptibility, 
which allows the measurement of local penetration depths, spins, and superconducting 
fluctuations with unprecedented sensitivity and resolution; MFM imaging of Meissner 
force derivatives, which allow spatially resolved absolute measurements of penetration 
depths; and controlled manipulation of superconducting vortices. I expect to see many 
other exciting developments in the future. 
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